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A. STATEMENT OF PROBLEM 

As microelectronic components continue to decrease in 
size and increase in component density, the problem of heat 
dissipation becomes increasingly more critical in the design 
of circuits and components. A five-order-of-magnitude 
increase in circuit integration in the last 25 years has 
been associated with successive revolutions in components 
technology. 

With the recent development of CMOS, industry has seen 
an almost three-order-of-magnitude reduction in size, an 
order-of-magnitude increase in the characteristic component 
dimensions and perhaps, most importantly, a precipitous drop 
in transition switching energy from 10-7 J in 1960 to nearly 
10713 3 in present day devices. It is anticipated that 
power dissipation of 100 W in a 10 million component chip 
could be encountered by 1990. This would result in 
component heat fluxes in excess of 100 W/cm? [(Bar-Cohen, 
nice pp. 159-175). The strong requirement to keep 
components at low temperature levels, 100°C to 125°C, does 
exist, since for every 20°C decrease in components 
temperature, the components failure rates are cut in half 


fOktay, 1986, pp. 36-42]. 


B. BACKGROUND 

Air cooling is the standard design in today's systems. 
Advances in air cooling technology now include impingement 
cooling and combined conduction and air cooling technology 
which offers even higher cooling rates than natural or 
forced convection air cooling. However, air as a heat 
transfer medium is extremely limited not only in its ability 
to transfer heat but also in its capacity to carry the 
transferred energy with a reasonable temperature increase 
(Baker, 1972, DPpeaee l3=22293 

Special component surfaces have been investigated. 
Features such as Turbulence Promoting Fins (TPF) have 
Significantly higher cooling enhancements when used with 
special gases such as Helium as the closed cycle coolant 
(Kishimoto et al., 1983, pp. 286-93]. By doing so, a 300% 
greater dissipation rate can be reached than in an open 
cycle with air as the coolant. 

Forced convection in liquids offers greater heat 
transfer coefficients than does natural convection in 
liquids or air. However, forced convection introduces two 
major problems: (i) greater amounts of expensive coolant 
are required and (ii) pumps and external plumbing increase 
the size and reduce the reliability by introducing the 
possibility of vibration related failures to the system 


(Baker, 1973, pp.) 63-1725 


Bolling can produce very large scale heat transfer 
coefficients, but it also introduces major problems: (i) 
boiling restricts the physical design of the system, (ii) if 
the surface is very smooth, as most electronic components 
are, a high degree of superheat may be required so that 
mucleation can occur. The occurrence of periodic boiling 
due to the high local heat flux may result in severe local 
turbulence and mechanical damage and (iii) boiling may 
result in the concentration of impurities at the nucleation 
site (Baker, 1973, pp. 163-173}. 

C. STUDIES ON NATURAL CONVECTION COOLING OF ELECTRONIC 

DEVICES 

1. Numerical Studies 

Of the different heat transfer modes mentioned 
above, immersion cooling by natural convection offers not 
only high heat transfer potential, but also a practical 
means of controlling circuit temperatures. However, only 
limited heat transfer data are available for immersion 
cooling of discrete, protruding heat sources [Liu et al., 
1987, pp. 309-30]. Early studies of natural convection heat 
transfer were motivated by applications such as energy 
conservation and solar energy collection. The pertinent 
length scales were several order of magnitude higher then 
those applicable in cooling of electronic equipment. 


However, recent studies have focused on computational and 


experimental studies of natural convection cooling in 
electronic components. 

(Han et al., 1984, pp. 21-28] have used numerical 
simulation to investigate the transient laminar natural 
convection in a rectangular enclosure with a localized 
constant heat flux source. [Acharya et al., 1983, pp. S335 
66] have investigated the natural convection phenomenon in 
an inclined square box shaped enclosure which has energy 
sources in it. ([Jaluria, 1985, pp. 883-92] investigated the 
interacting natural convection flows generated by flush, 
isolated thermal energy sources, such as electronic 
components, located on a vertical adiabatic surface. 
[Jaluria, 1982, pp. 223-27] also analytically studied the 
natural convection flow due to multiple isolated heated 
elements located on a vertical surface. (Torok, 1984 73ipp- 
49-57] used the finite element technique to predict 
component and board temperatures and heat fluxes. 

(Kuhn et al., 1986, pp. 55-62] have investigated the 
transient three-dimensional natural convective flow and heat 
transfer within a rectangular box following a sudden 
temperature change on a rectangular portion of one vertical 
wall. They found that when two elements are symmetrically 
placed on the wall, the distance between them has only a 
slight effect on the mean heat transfer distribution. 
However, it does have a significant effect on the local 


Nusselt number. (Shakerin et al., 1986, pp. 17-22] 


investigated the natural convection flow adjacent to a 
heated wall with single and repeated, two-dimensional, 
rectangular roughness elements. fOosthuizen et al., 1987, 
pp. 1-7] studied the heat transfer rates from an element of 
Square cross-sectional shape held at a uniform temperature 
and mounted on a vertical adiabatic wall of a tilted square 
enclosure. They found that the mean heat transfer rate from 
the block did not vary continuously with the angle of 
inclination. 

(Liu et al., 1987, pp. 486-97] have conducted a 
finite difference numerical study of natural convection flow 
in a rectangular enclosure of an array of simulated 
components that are mounted on a vertical wall. The results 
of this numerical study have been presented for enclosure 
widths of 18 and 30 mm. It was concluded that: 

- The temperature field in the enclosure is characterized 
by boundary layer regions surrounding individual 
components. 


- There are only small interferences among the components, 
especially for the lower rows. 


- The maximum temperatures on the components surfaces are 
located on the upper horizontal surfaces of the 


components. Lower temperatures have been found on the 
lower horizontal faces of the components in the bottom 
row. 


- The results for the 30 mm width are essentially the same 
as those for the 18 mm, indicating only a minor 
dependence on the enclosure width. 


fuiu et al., 1987, pp. 309-30] have also observed a local 


oscillatory surface temperature response in the _ same 


enclosure. Their study found that the components surface 
temperatures have a tendency to fluctuate within a range of 
approximately 3°C with a period of 4 seconds. 

2. Experimental Studies 

Numerous experimental studies have been conducted 
concurrently with the numerical studies described above. 
Two of the earliest studies were conducted by [Baker, 1972, 
pp. 213-22] and (Baker, 1972, po wic 7 The purpose of 
these investigations was to examine liquid immersion 
techniques for cooling miniature heat sources, specifically 
effects of component size on the convective heat transfer. 
[Goel et al., 1986, pp. 1341-46] experimentally investigated 
the two-dimensional natural convection flow from a long, 
finite size heat source located on a vertical or inclined 
plate in air. The principal items of interest were the 
nonboundary layer effects that arise in the flow, 
particularly in the vicinity of the source. Flow separation 
was observed at large inclinations. 

([Milanez et al., 1986, pp. 1347-52] experimentally 
studied the interaction of two horizontal, wide strips ona 
vertical adiabatic wall in both air and water. They found 
the results of their experiment in good agreement with past 
numerical studies. [Ortega et al., 1986, pp. 5-15) 
presented the results for a series of experiments in air 
designed to measure the natural convective heat transfer 


rates and the thermal characteristics of a rectangular array 


of cubical roughness elements mounted on an insulated plane 
wall, with and without a shrouding wall. [Parks et al., 
1987, pp. 91-109] simulated microelectronic circuits with 
thin foil heaters for both flush and protruding components. 
They found that the protruding heaters have approximately 
15% higher coefficients and that the upper heaters have 
varying coefficients depending upon the array configuration. 

A number of studies have been conducted at the Naval 
Postgraduate school. (Knock, 1983] studied the effect of 
the location of a single protruding heater in an enclosure 
filled with water. (Pamuk, 1987] investigated the heat 
transfer characteristics of an array of nine protruding 
components immersed in a dielectric fluid in an uninsulated 
chamber. (Hazard, 1987] studied the natural convection 
liquid cooling of a vertically arranged components array 


with and without a shrouding wall. 


D. OBJECTIVES 
This experimental work was motivated by the numerical 
study of [Liu et al., 1987, pp. 486-97] and the recent 
experimental studies of [Pamuk, 1987]. The objectives for 
this thesis are: 
- To design and build a circuit card of nine simulated 
components and determine the rate of heat transfer for 


various power inputs. 


~ To determine the natural convection flow patterns 
developed in the chamber by means of flow visualization. 


~ To correlate the data in terms of dimensionless heat 
transfer parameters such as Rayleigh vs. Nusselt number. 


IIT. EXPERIMENTAL APPARATUS 


A. GENERAL CONSIDERATIONS 

This experimental study is a continuation of past 
experiments conducted at the Naval Postgraduate School 
(Knock, 1983; Pamuk, 1987; Hazardyaeo2 ie The purpose of 
this work was to examine natural convection heat transfer 
and fluid flow from a 3 x 3 element array of simulated 
electronic components. A major change from the previous 
experiments was the design of the simulated circuit card. 
In [{Pamuk, 1987] an arrangement of vertical components was 
used. The present study utilized a horizontal arrangement. 
A decision to increase the card thickness was made to reduce 
the contamination of the dielectric fluid by the loose 
insulating material previously used and also to reduce the 
conduction loss through the card. 

The experimental apparatus consists of a chamber filled 
with a dielectric fluid. Top view of the chamber is 
provided in Figure 2.1. The top and bottom surfaces of the 
chamber are maintained at a prescribed temperature by means 
of two heat exchangers and two separate chilled water 
circulation baths. One vertical wall of the chamber 
Simulates a circuit card. It contains nine simulated 
components ina 3 xX 3 array, while all other vertical walls 


are insulated. The geometrical arrangement of the simulated 
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components on the card is seen in Figure 2.2. The 
components were chosen to simulate an array of 20 pin Dual- 
inline-packages (DIPS). Their dimensions are identical to 
those used by [Liu et al., 1987, pp. 486-97; Liu et al., 
1987, pp. 309-30; Pamuk, 1987] in order to permit future 
comparisons of experiments and computations. 

The heating elements were made of aluminum and heated by 
means of small foil type resistance heaters attached to the 
bottom surfaces. The temperatures of the faces of the 
elements were measured by means of embedded thermocouples. 
Data were recorded and processed uSing an automatic data 
acquisition system. The major components of the system are 


described in the following section. 


B. COMPONENTS 
1. Heating Element 
Each heating element was a block of aluminum 6mm by 
8mm by 24mm (Figure 2.3). The temperatures at the faces of 
the blocks were measured by means of thermocouples placed 
within grooves (0.5mm deep) on each surface. The front 
surface thermocouple was passed through the block by means 
of a 1.02mm hole drilled through the block. All grooves 
were filled with high thermal conductivity epoxy and 
smoothed prior to attachment to the simulated circuit card. 
2. Thermocouples 
All thermocouples were 10-mil (0.254mm) copper- 


constantan. These small wires were chosen to minimize the 
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material removed from the block. The thermocouples were 
manufactured by using an arc welder designed for such 
purposes. The setting for both current time andé gas 
pressure (argon) was kept at 'medium.' Due to the extremely 
small diameter of the wire the gas pressure was carefully 
regulated to less than 5 psig. The thermocouples were also 
checked for continuity due to the fact that the arc welder 
tended to embrittle the copper wire at the bead. 
3. Heaters 

A nearly constant heat flux condition was maintained 
at the base of the block by use of thin (0.18mm) foil type 
heaters (Figure 2.4). Each of the heaters was measured to 
have a resistance between 10.5 and 10.6 ohms. These heaters 
were designed to align with the blocks. Holes were provided 
within the heaters to pass the various thermocouple wires. 

The power leads of the heaters were gold plated to 
ensure a low electrical contact resistance. Wires from the 
power supply were attached to the leads by means of 
soldering. After ensuring proper alignment, the heaters 
were bonded to the blocks by means of a high thermal 
conductivity epoxy. 

4. Simulated Circuit Card 

Three grooves of 0.5mm depth were made on the front 
of the card in order to allow the heaters to be mounted 
below the surface and thus prevent the dielectric fluid from 


seeping under the blocks and introducing a variable, unknown 
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resistance. The heating elements were bonded to the circuit 
card with a low thermal conductivity epoxy. 
5. Test Chamber 

The test cnamber was constructed of plexiglass (13mm 
thick) and had inside dimensions of 120mm in height, 144mm 
in length, and 30mm in width. The wires were fed through 
the grooves in the back of the circuit card to a "Tygon" 
tube on the back wall. The chamber also contained two small 
pressure equalization outlets, located at the top and the 
bottom of the front vertical surface. 

6. Heat Exchangers 

The two aluminum heat exchangers measure 38mm by 
65mm by 274mm. The details of the heat exchangers are 
avallable in [Knock, 1983]. Aluminum plates 3mm thick were 
used as top and bottom surfaces. These plates were grooved 
2.5mm to provide slots for the thermocouples. This left a 
0.5mm wall thickness of aluminum between the thermocouple 
bead and the surface in contact with the dielectric fluid. 
The resulting conductive thermal resistance was found to be 
negligible. The heat exchangers were each attached to 
individual coolant baths by means of "Tygon" tubing. Four 
10 mil (0.254mm) thermocouples, two on the top and two on 
the bottom, were used to read the heat exchanger 


temperatures. 
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7. Flow Visualization 

Visualization was accomplished with magnesium 
powder, 325 mesh. Magnesium has a density of 1.92 g/cm? 
compared to 1.76 g/em? for FC-75, resulting in a large 
suspension time of the particles. The beams from the light 
sources were spread into plane sheets which illuminated the 
particles suspended in the FC-75. This technique allowed 
for the visualization of a Single plane of the flow field. 
Other planes were examined by the realignment of the 
apparatus. A four milli-watt helium neon laser and a 
cylindrical lens were used for the time lapse camera photos 
(Figure 2.5). 

8. Assembly 

After sliding the card into the slot in the chamber, 
both the thermocouple and power wires were run out through 
the 1.27cm diameter extension pipe on the back of the 
chamber to the connected "Tygon" tubing of approximately 1.0 
meter in length. 

Rubber o-rings were placed in grooves located on the 
top and bottom surfaces of the chamber and the thin aluminum 
plates were set in place. The heat exchangers were then set 
on the plates and clamped by means of specially made 
aluminum clamps to ensure a proper seal. The inlet and 
outlet tygon tubes were attached to the constant temperature 


baths for each heat exchanger. Finally, the pressure 
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schematic of the entire assembly is shown in Figure 2.6. 
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Tit. EXPERIMENTALSCER CC ERORE 


A. APPARATUS PREPARATION 

After assembling the apparatus, the chamber was filled 
with the dielectric fluid through the top vent tube by means 
of a small funnel and a hypodermic syringe. To damp 
disturbances, the entire apparatus was placed on a separate 
platform. This platform also served as a means of leveling 
the chamber by use of adjustable screw legs. Since the 
entire apparatus was sealed by means of four clamps it could 
be tilted in order to remove the last remaining air bubbles 
inside the chamber. 

The use of a single constant temperature bath providing 
cooling water in series to the top and bottom heat 
exchangers proved inefficient to maintain both the top and 
the bottom heat exchangers at the same constant temperature. 
A second bath was added so that each heat exchanger was 
supplied individually and could be adjusted separately. The 
ten terminals on the power supply panel were used to supply 
power to the heaters on the simulated components. The first 
terminal was the main input, the next nine were the 
terminals to the individual heaters located on each 
component. These heaters were connected in series with the 


2 ohm precision resistors. This arrangement was repeated 


20 


for each component. These nine units were then connected in 


parallel to the main input. 


B. DATA ACQUISITION 

Prior to energizing the power supply or the ponetane 
temperature baths all the thermocouples were scanned to 
measure the starting steady state temperatures. The 
constant temperature baths were set at 10°C by means of 
adjustment dials located on the baths. Throughout each 
experiment the bath supplying the top heat exchanger 
required constant monitoring. As the power levels 
increased, the bath temperature was reduced accordingly to 
maintain the top heat exchanger thermocouple readings at a 
Constant 10°C. 

The desired voltage was obtained by setting the power 
supply to the desired level. The reading on the power 
supply was verified by means of reading channel number 61, 
the main power input channel, on the data acquisition panel. 

Approximately four hours were allotted for the system to 
achieve steady state. For this experiment, the steady state 
condition was defined as a majority of the thermocouple 
readings changing less than 0.2°C during three separate 
readings taken at ten minute intervals. The bath 
temperatures were checked to ensure they were at 10-@e If 
they were in agreement then the data was recorded by means 
of the program "ACQUIRE" (Appendix C). This program scanned 


all 76 channels and converted the emf signals from the 


Pa. 


thermocouples to temperatures in degrees Celsius. The power 
to the heater was determined by utilizing the precision 2 
ohm resistor. This resistor was connected in series with 
the heater and allowed the measurement of the current 
through this configuration. This current was’7 then 
multiplied by the measured voltage of the heater to produce 
the power input to the component. 

A second program "Fastscan" (Appendix C) was developed 
to scan three individual channels. Each of the three 
channels, numbers 1, 13, 32, was scanned individually at a 
rate of 3 scans per second for a period of 200 seconds. The 
main intention in doing so was to isolate any periodic 
behavior caused by the flow past these three thermocouples. 
A third program "Plot" (Appendix C) was developed to plot 
these results on a HP model 7470A plotter. 

Flow visualization was conducted by means of 325 mesh 
Magnesium particles and a four milliwatt laser focused with 
a cylindrical lens. The magnesium particles were injected 
into the chamber by way of the bottom vent tube. A motor 
driven Nikon F3 camera with a Micro-Nikkor 55mm lens and a 
Nikon MT-2 Intervalometer was used for the long-time 
photography. The time exposure photographs of length 8 or 
20 seconds were taken at the various power levels where heat 
transfer data had been obtained. The photographs were taken 
of the three columns of components by placing the camera at 


a ninety degree angle to the plane created by the laser. 


Ze 


This configuration was then reversed to obtain photographs 


of three separate planes in the chamber. 


oF DATA ANALYSIS 
1. Determination of the Nusselt Number 


The Nusselt Number is defined as the following: 


where L is a characteristic length, taken here as the length 
of the component in the vertical direction, kg is the 
thermal conductivity of the fluid, and h is the convective 
heat transfer coefficient. All fluid properties were taken 
to be a function of the average film temperature. 

To evaluate the heat loss, one thermocouple was 
placed directly behind each of the heated blocks. A series 
of nine thermocouples were also attached to the back of the 
Sage« These temperatures were used to calculate the 
conduction loss through the back of the card. 

One dimensional conduction resistance, Raon, for the 


simulated board was determined from: 


3 
Roon = ka (Si2) 


where L is the length along the conduction path i.e., the 
thickness of the board, k is the thermal conductivity of the 


material (Appendix A) and A is the area normal to the 
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conduction path. The area normal to the path of conduction 
was calculated as the area directly behind the components 
back, 192mm°. By doing so spreading losses were neglected. 


The conduction loss is then: 


pi - T, 
Qcon = R (3.3) 
con 
Power input to each component is calculated from: 
(Va. = Vay 
un) ho ob 
Qin = a (345) 


where Vj, is the input voltage applied to both the 2 ohm 
resistor and the heater which is in series with it. V} is 
the voltage drop across the heater. The difference between 
Vin and Vp, gives the voltage drop across the 2 ohm resistor. 

The net energy added to the fluid in the chamber is 
the difference between Qi, and Quon which by Newton's Law of 


cooling iss 


Qconv = Qnet = HA(Tayg - To) 


Qin ~ Qcon (3.5) 


where, 


The total wetted surface area of the 
components, 


> 
il 


Tavg = Area weighted average surface temperature, 
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calculated by, 


iL 
Tavg = “aA (3.6) 


where, 
T; = Temperature of the ith surface of the 
components, 
A; = Area of the it? surface of the components. 
The summation of the index i varies from 
Tce 5:; 
helen 
To = Average heat exchanger temperature. 
Thus, 
Q 
oa oa oes 
EOG wavG G 


The above procedure was repeated for each of the 
nine components at each power level. The Nusselt Number was 


then calculated from Equation 3.1. 


2. Determination of the Rayleigh Numbers 
The definition of the temperature based Rayleigh 


number is: 


Rat = Gre * Pr (a7) 
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where the Grashof Number (Gry) is defined as, 


3} 
Silos (AE ed 
Grp = —__-“9 = (3.8) 
v 
and the Prandtl Number (Pr) is defined as, 
anil 

Pree (32) 
where, 

5 ef e 

See cue) 
" Pp 


All the properties in the above equations are 


calculated at the average film temperature: 


a + To 
Tr = g 5 C3 Op) 
Therefore: 
gBL(T, - 7) 
Rat = ae ( 3..0518) 


3. Determination of the Flux Based Rayleigh Number 


The definition of the flux based Rayleigh number is: 


Rag = Gre * Pr (3°. W223) 
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all film properties are calculated at the average film 


temperature as before and: 


4 
gBL Chet 


Crp = aS Ls) 
Sau Atot 


Functional relationships for the properties at the 
film temperature are given in Appendix A, together with 


sample calculations. Uncertainty calculations are given in 


Appendix B. 
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IV. RESULTS 


A. FLOW VISUALIZATION 
Flow visualization was accomplished using Magnesium 
particles as a flow marker and a Helium Neon laser to 
illuminate the particles. The results provided insights 
into the flow patterns existing at various power dissipation 
levels. As expected, the particle traces suggest highly 
three-dimensional flow patterns at all power levels. 
Particles can be seen entering and leaving the plane of 
light in almost all the long time exposure photographs. A 
multitude of flow patterns was observed for the wide range 
of input power levels in the of 0.1 to 3.1 watts. 
Representative flow visualizations are presented in Figures 
4.1 to 4.14. 
1. Flow Patterns at Component Midplanes 
The following flow patterns are in a vertical plane 
through the center of components 4, 5, 6 (Fig. 2.2). The 
laser was placed in front of the chamber and the camera 
located to the side (Fig. 2.5). The camera exposure time 
was set at 10 seconds. The three dark regions located on 
the left side of the photographs are the three rows of 
components as seen in side view. The top and bottom of the 
chamber and the horizontal walls can also be seen in the 


photos. 
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Figure 4.1 is the flow pattern for a power input of 
0.1 watts and indicates a pattern that contains two 
individual eddies located between the rows of components 
along with an overall circular pattern that encompasses the 
entire chamber. A boundary layer like flow is seen in the 
vicinity of the heated components. A dark region located at 
the bottom of the photograph is evident in all photos taken 
throughout the experiment. This region is the stagnant, 
dense cold fluid layer in the vicinity of the bottom heat 
exchanger. This region proved to be impermeable at all 
power levels. 

Figure 4.2 for power input of 0.3 watts indicates 
the same overall pattern with the exception of a 
constriction of the overall circular pattern at the 
approximate height of the second row of components. The 
slight increase in the upper eddy size is also evident. 

Figure 4.3 for a power input of 1.5 watts indicates 
an increase in the thickness of the boundary layer emanating 
from the second row of components along with increased 
entrainment of the particles located below the second row of 
components. Also evident is the disappearance of the two 
distinct eddies seen earlier and strong evidence of three 
dimensionality in the flow as seen by the crossing of the 


particle paths. 
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Power @ 0.3 W 
Side view 


Figure 4.2 


Power @ 0.1 W 
Side View 
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Figure 4 
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Power @ 1.5 W 
Side View 


BUEguire:. 4. 3 


oat: 


Figure 4.4, for a power input of 2.7 watts, shows a 
buoyant fluid layer adjacent to the components and random 
motion throughout the remaining chamber. Also seen are 
entrainment of particles at the approximate height of the 
second row of components, and evidence of aan area, 
immediately above the black region at the photo's botton, 
where particles are basically suspended. 

Figure 4.5, at 3.1 watts, demonstrated a downward 
flow pattern in the vicinity of the bottom row. This 
downward flow was forced past the face of the bottom chip in 
this® eelumn-. The flow also demonstrated a perpendicular 
entrainment of particles between the first and second row of 
components, and a large increase in the boundary layer at 
the third row of components. 

Figures 4.6 and 4.7 show flow patterns in the 
vertical plane at the mid-span of components 7, 8 and 9. 
These two photographs were taken approximately 5 seconds 
apart at a power setting of 3.1 watts. They demonstrate the 
strongly time-dependent behavior of the flow. In Figure 
4.6, entrainment into the upward moving buoyant layer is 
seen to be present almost as far down as the bottom of the 
lowest component. Far away from the components the core 
fluid exhibits no preferential direction of motion. In 
Figure 4.7 the flow pattern has changed. Near the lowest 
row of components the flow now appears almost quiescent. 


Also, the buoyant region appears thinner than in Figure 4.6. 
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Figure 4.4 Power @ 2.7 W Figure 4.5 Power @ 3.1 W 
Side View Side View 
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Figure 4.6 
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2. Flow Patterns in Other Planes 

Figure 4.8 is a schematic of the three planes that 
were investigated when the positions of the camera and the 
laser were reversed. Figures 4.9, 4.10 and 4.11 are at a 
power setting of 0.1 watts and demonstrate the flow patterns 
in the back, middle and front planes, respectively. 18 e 
Figure 4.9, the back plane, the flow appears to be somewhat 
like a large circular pattern that is bounded by the 
physical dimensions of the chamber and exhibits no 
noticeable flow directly up the center of the chamber. 
Figure 4.10, the middle plane, indicates a strong upward 
flow in the center of the chamber and two areas of eddy 
formation located approximately mid-height in the chamber. 
Finally, Figure 4.11, the front plane, demonstrated 
approximately the same behavior as the middle plane; a 
general area of motion in the center of the chamber along 


with eddy formation areas on either Side. 


' SANT Id 





Figure 4.8 Top View of Planes 
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Figure 4.10 Middle Plane @ 0.1 Watts 





Figure 4.11 Front Plane @ 0.1 Watts 





Figures 4.12, 4.13 and 4.14 are the same three 
planes at a power setting of 3.1 watts. Figure 4.12, the 
back plane, indicates regions of upward flow near the 
components. These patterns characterize the boundary layer 
nature of transport in the vicinity of the elements. 
Relatively short particle traces in regions away from the 
components appear to result due to the strongly three- 
dimensional transport. The dark area in the very bottom of 
the photograph is again pointed out to demonstrate the cold 
isothermal layer developed by the bottom heat exchanger. 

Figures 4.13 and 4.14 are views of the center and 
front planes respectively. The flow in both planes appears 
strongly time dependent. Figure 4.14 shows relatively 
longer particle traces than Figure 4.13, indicating stronger 
three dimensional effects near the components. 

These 13 photographs illustrate the complex three- 
dimensional, time dependent flow patterns found in this 
study. Obtaining a detailed understanding of these effects 
is essential in achieving predictive capability for the heat 


transfer in these configurations. 


B. REDUCTION OF HEAT TRANSFER DATA 
The temperature and reduced data tables corresponding to 
the eight power levels are as follows: (see Figure 2.2 for 


identification of the component numbers). 
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Figure 4.12 Back Plane @ 3.1 Watts 
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Figure 4.13 Middle Plane @ 3.1 Watts 
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Figure 4.14 Front Plane @ 3.1 Watts 
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TABLE 1 


TEMPERATURE DATA FOR INPUT POWER Oiwe= Oo WW 


Teese NESULI3S ARE =fORED IN FILE: 98098133 
AME SET TEMP LAS: 2k 
VOLIMETER SEADING wAS: 2.9832 et! 
BATH TEMP 46S: | Oe 


THIS RUN WAS RECORDED WITH FLOW VIZ 
eV TEMPERATURES ARE IN DEGREES CELCIUS 


CENTER TOP RIGHT GEE T corenoM 

CHIP NO!: 16.9546+00 16.336+90 16.15F+00 16.156+00 !16.118+90 
POWER (WATTS): 29.9325-02 

CHIP NO2: 16.686+00 1!6.31FE+00 16.346€+00 16.19&*00 16.30&+00 
POWER (WATTS): 29.94E-02 

CHIP NO3: 16.74€*00 16.226+00 !6.466€+00 16.33E6+00 16.636 +00 
POWER (WATTS) Fa S| 8).5\), 10 

CHIP NO4&: 16.26€+00 16.23£+00 16.086+00 %15.82E+90 16.02E6+00 
POWER (mATTS) (ZS) SVS) a Ut 

CHIP “0S: 16.526€+00 16.:S&+d0 %16.!12£+00 i6.056+00 16.306+00 
POWER (WATTS): 29.672-02 

CHIP NO6: 16.51€+00 16.50€E+00 16.33£+G0 16.33£+00 16.58E&+*00 
POWER (WATTS): 29.€64E-02 

CHIP NO7: 16.646E*00 16.42£+00 16.12£+00 16.12E+00 16.32E+00 
POWER (WATTS): 30.0!2-02 

CHIP NOS: 16.716€+00 16.42£+00 16.40€*00 16.40€+00 16.28E€+00 
POWER (WATTS): 30.12£-02 

CHIP NOS: 16.55£+00 16.286€+00 16.236+00 16.36£+00 16.41£+00 
POWER (WATTS 29.37E-02 
HEAT EXCRANGERS TEMPERATURES: RIGHT LEFT 

BGrron; ° 99.12£-91 99,39E-01 
TGR? 10.386+00 938.935E-0! 

BACK PLANE TEMPERATURES ARE : 


T(SS): 14.67&+00 
T(S6): 15.41E&+00 
T(S7): 15.276 +00 
(72): 15.S7E+00 
™(73): 15.556 +00 
T(74): 16.06E6+00 © 
(75): 15.56c+00 
TO76):) «615.5626 +00 
T(77): «615.456 +00 
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BACK 


fe J5c © OU 


ire > OE. +00 


17.40E +00 


We ee re) ee O19, 


emer su 


18.156 +90 


Woe tO 


18.02E+00 


17.635E+00 


FeokE 2 


TEMPERATURE DATA FOR INPUT POWER Qj, = 0.7 W 


CHIP "I01: 
POWER 
CHIP NO2: 
POWER 


CHIP NO3: 
POWER 


CHIP NO4: 
POWER 


CHIP NOS: 
POWER 


CHIP NOG: 
POWE 
CHIP NO7: 
POWER 
CHIP NO8: 
POWE 


GHre NOS: 
POWER 


AMBIENT TEMP WAS: 
VOLIMETER READING HAS: 


BATH TEMP WAS: 
THIS RUN WAS RECORDED WITH FLOW VIZ 


THESE RESELIS ARE GIURee meee. 


C37 


Lose 


Dgdvxe2 10 


Oe oie We 


ALL TEMPERATURES ARE IN DEGREES CELCIUS 


CENTER 


25156 <0 
(WATTS): 


23.266 +00 
(WATTS): 


c3.16E+00 
(WATTS): 


2euc se “00 
(WATTS): 


23.11E+90 
(WATTS): 


22.69 +00 


R (WATTS): 


23.85E 400 
(WATTS): 


23732700 


R (WATTS): 


23.07£ +00 
(WATTS): 


HEAT EXCHANGERS 
BOTTOM: 


BACK 


TOP: 


TOP | 


22550. +00 
71.08&-02 


22.3/&+00 
Maa S.3)74 


coe 0 
FUSE sue 


22.426 +00 
/0.18E-02 


FeV Esha aol) [) 
70.47E-02 


22.77E+00 
70.38E-02 


22.76E+00 
Wecoe 02 


22.63E +00 
11 e5se 202 


22.54E+00 
70.352 -02 


RIGHT 
PafeISVE AN LD 


22.43E +00 
22755e700 
22-S1Es00 
21.94E+00 
22.!186+00 
D2 E200 


22.47E +00 


22 -05e 700 


etal 
22. 396000 


22.09&+00 
22. 14E +00 
21.69E€+00 
21.85E+00 
Pee E00 
222 15e +00 
22, Soe uU 


A EMI 


BOTTOM 


22 


TEMPERATURES: 


PLANE TEMPERATURES ARE : 


WSs = 
T(S6): 
WSs 
Wz2 ye 
mG 30: 
T(74): 
T1759): 
T7694 
l(7e 


es at “8 oH OS 8 SOD ot =O 
MDD DMWW ON 


".67E+00 
. 766 +00 
.0SE+00 
.08E +00 
.90E+00 
116 +00 
.02£+00 
.79E +00 
.71£+00 


“a 


RICHT 
Sele slat) ) 
3e ce Ot 
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Ee | 
98.46E-01 
33.598E-01 


~24& +00 


soe 00 


-73E +00 


. 166 +00 


27E +00 


.60£ +00 


-39E +00 


~47E +00 


-20E +00 


24 


26 


24 


26 


25 


BACK 


.06E+00 


.05E +00 


-58E +00 


- 16£ +00 


7 toe +00 


.1Ze+0u 


2554200 


- (7E+0G 


~41E +00 


ieee 3 


TEMPERATURE DATA FOR INPUT POWER Con see ~6W 


TRESe RESULTS ARE STORED IN FILE: 080990745 
BGG Tec! : TEMP Was: eo AG 
VOLIMETER FEADTNG WAS: 4.1107 bl 
Bate TEMP ane: (ee 


THIS RUN WAS RECORDED WITH FLOW VIZ 


ALL TEMPERATURES ARE IN DEGREES CELCIUS 


CENTER TOP RIGHT LEFT BOTTOM 
CHIP NOt: 28.626 +00 27.2356+N0 27.336+00 27.16€+00 26.92£+00 
POWER (WATTS): 11.S57€-01 
CHIP NO2: 28.73£+00 27.05€+*00 27.13£+00 26.72E+00 27.11£+00 
POWER (WATTS): 11.S8£-01 ; 
CHIP NO3: 28.47£+00 27.17£+00 27.49£+00 26.82E£+00 27.69£+00 
POWER (WATTS): 11.S4E-01 
CHIP NO3: 28.47£+90 27.046€+90 27.146£+00 26.138E+00 26.95E6+00 
POWER (WATTS): 11.438-0! 
CHIP NOS: 28.21£+00 26.29&+90 26.42£+00 26.30E+00 27.03£+00 
POWER (WKATTS): 11.48€-01 
CHIP NOG: 27.45SE+00 27.90£+00 26.98E+00 27.00£+00 27.76&+00 
POWER (WATTS): 11.46€-01 
CHIP NO7: 29.12E+00 27.77£+00 27.03£+00 26.97E+00 27.25£+00 
POWOR (WATTS): 11.605£-01 
CHIP NOS: 28.S6£+00 27.4SE+00 27.27£+00 27.40£+00 27.18£+00 
POWER (WATTS): 11.65c-01 
CHIP NO9: 28.13£+00 27.32£+00 26.60£+00 27.12E£+00 26.72E£+00 
POWER (hATTS): 11.55e-01 
HEAT EXCHANGERS TEMPERATURES: RIGHT Bean 
BOTTOM: . 99.44£-01 99.77E£-01 
TOP: 10.05e*00 936.50E-0! 
BACK PLANE TEMPERATURES ARE : 
1055): 19.74E +00 
T(S56)3 21.13E6+00 
T(S7): 21.79&+00 
T(72): 21.42E£+00 
T(73):3) 21.346 +00 
T(74): 21.418 +00 » 
Waste a2 le Sce un 
T(76): 20.93£+00 
77 se coe ocee 700 
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BACK 
eoeove “00 


S17438+00 


30.65£ +00 


SOF Ceo 


30.00£+00 


33.¢8E&+00 


30.656 +00 


Sanc UE +00 


Sane +00 


TABLE 4 


TEMPERATURE DATA FOR INPUT POWER Qin = 1-5 W 


CHIP NOI: 
POWE 
CHIP NO2: 
POWE 
CHIP NO3: 
POWE 


CHIP NO4: 
POWER 


CHIP NOS: 
POWER 


CHIP NO6: 
POWE 


CHIP NOQ?7: 
POWER 


CHIP NO8: 
POWE 


CHIP NOS: 
POWER 


THESE RESULTS BRE = SI CREo erie ee. 


AMBIENT TEMP WAS: 
VULIMETER FEADING wAS: 
US 


BATH TEMP «AS: 
THIS RUN WAS RECORDED WITH FLOW VIZ 


ae ale 


4.5442 


0803992015 


eh} 


ALL TEMPERATURES ARE IN DEGREES CELCIUS 


CENTER 
S201 +00 


R (WATTS): 


ge. 2ce 200 


R (WATTS): 


oa. 366 700 


R (WATTS): 


32.60E +00 
(WATTS): 


32.28E£+00 
(WATTS): 


30.86E +00 


R (WATTS): 


33.43E+00 
(WATTS): 


32.822 700 


R (WATTS): 


Seecre uo 
(WATTS): 


HEAT EXCHANGERS 
BOT : 


TOM: 
TOP: 


he 


Site Je +00 
149.77&-01 


30236 +00 
14,77E-01 


302252290 
14,72E-01 


30.86E+00 
14.58€-01 


Male flea )N 0) 
14.64E-01 


3l7 (e400 
14.02E~01 


31.83E+90 
14.80E-01 


31.45E+00 
14.87E-01 


Sre0ce +00 
14.74E-01 


TEMPERATURES: 


RIGHT 
g12C3e +00 


Si 02E «00 


Bh Fd Se 20) 


31.04E+00 


30.12&+00 


30.66&+00 


30.89E+00 


Ste e7eod 


30.36E+00 


RIGHT 
92. 49E=07 


CEG 
31.088 +90 


30.42E +00 
30.50E +00 
29). 86°00 
295228400 
30.72E+00 
30.81&+00 
31.40&+00 


30.37E+00 


BERT 
eo sOdaeut 


BOTTOM 
3037 8e+00 


312076200 


31.66e+00 


3077736700 


JO, ene +00 


31.43E+00 


3IVe1Z2E+00 


S31. 226400 


30.306 +00 


22.0 Ie-0 ego. SoU. 


BACK PLANE TEMPERATURES ARE : 


T($5): 92! 
T(56): 23 
T(S7): 24, 
EC72)3- “Zee 
T(73)2: 2 
T(74): aye 
1075) 2am 
T(76): 2 
77) ee 


~87E +00 
~45E +00 


S36 +90 
B4E+00 


-52E+00 


83£+00 | 
10£+00 


-cSE+00 
~ 10E+00 
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BACK 
34.47E£ +00 


36 .44E+00 


35-590 700 


34.86E +00 


34.66E*00 


39.92E+00 


35.45£+00 


38.77E+00 


372258200 


TABLE 5 


TEMPERATURE DATA FOR INPUT POWER Qin =1.9W 
THESE RESULTS ARE STORED IN FILE: ODe8t001139 
AMBIENT TEMP WAS: PLR 
VOLTMETER READING WAS: $.2644 eh] 
EATH TEMP WAS: Las 
THIS RUN WAS RECORDED WITH FLOW VIZ 
ALL TEMPERATURES ARE IN DEGREES CELCIUS 
CENTER TOP RIGHT Ceri BOTTOM BACK 
CHIP "i011: 37.4SE+N00 35.42E*0N 35.24€+00 35.29F+00 34.90E+00 39.63EF +00 
POWER (WATTS): 19.10E-03 
CHIP NO2: 37.466£+00 34.63E*00 35.006+00 34.29E+N0 34.99£&+00 42.43E£+00 
POWER (WATTS): 19.131E-013 
CHIP NO3: 36.60€+00 34.89&6+N0 35.33E6+00 34.33E+00 35.S516+00 40.466 +00 
POWER (WATTS): 19.0SE-01 
CHIP NO4: 37.°22£+00 34.986+00 35.28FE+00 33.81E+*00 35.01E+90 40.20£ +00 
POWER (WATTS): 18.86E6-013 
CHIP NOS: 36.7NE*00 33.55€*N0 34.01&+*00 33.866+00 34.97E+00 39.836+00 
POWER (WATTS): 18.94E-01 
CHIP NO6: 35.0SE+*00 36.SS5E+00 34.S59E+00 34.9S5E*00 35.73£+00 48.38E£+00 
POWER (WATTS): 18.90E6-01 
CHIP NO7: 38.33E+00 36.17E+00 35.136+00 34.99£+00 35.34&+00 40.79E£+00 
POWER (WATTS): 19.1SE-01 : 
CHIP NOS: 37.38E+00 35.81E+00 35.42E+00 35.S7£+00 35.43£+00 45 .03£+00 
POWER (WATTS): 19.23E-01 
CHIP NOS: 36.47E+00 35.22E+00 34.366+00 35.106+00 34.40E+00 43.136+00 
POWER (WATTS): 19.07&-01 
HEAT EXCHANGERS TEMPERATURES: RIGHT seu 
BOTTOM: 99.34—E-01 10.026+00 
TOP: 95. .c9E-0! 92.47E-01 


BACK 


PLANE TEMPERATURES ARE : 


PSS): 23 
T(So): 


TUS7 eS. 
-81E+00 
25.38E +00 
~72E +00 
--9E+00 -, 


Troe: 
1(73): 
T(74): 
W7S ee te 


OLS ES Pace 
.85E +00 


T(77): 0 24 


»59E +00 
2 ooe +0 


S8E+00 


04E +00 
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TABLE 6 


TEMPERATURE DATA FOR INPUT POWER Qin = 2.3 W 


THESE RESULTS ARE STORED IN FILE: D81001740 
AMBIENT TEMP WAS: Cane 
VOLTMET&ER READING AS: 5§.8093 «#61 


BATH TEMP WAS: 1O2Ee 
THIS RUN WAS RECORDED WITH FLOW VIZ 


CHIP NO!: 
POWE 
CHIP NO2: 
POWER 
CHIP NQ3: 
POWE 
CHIP NO4: 
POWER 
CHIP NOS: 
POWE 
CHIP NO6: 
POWE 
CHIP NO?7: 
POWE 
CHIP NOB: 
POWE 


CHIP NOY: 
POWER 


ALL TEMPERATURES ARE 


CENTER 
42.48E +90 


R (WATTS): 


42.44E+1)0 
(WATTS): 


41.90E+00 


R (WATTS): 


42.09E +00 
(WATTS): 


41.49&+00 


R (WATTS): 


39.67E+00 


R (WATTS): 


43.53£ +00 


R (WATTS): 


42.48E+00 


R (WATTS): 


41.70E+00 
(HATTS): 


TOP 


39.94E +00 
23-07E-0; 


SJeeoe 20 
2s 06E 20" 


39.50E +00 
23, 00E-0) 


39.246 +90 
22./78E-01 


37.80E+00 
22.88E-01 


41.33€+00 
22.82E-01 


40.80E+00 
235126-01 


40.36E +00 
236220501 


39.86E+00 
23, 0235-01 


RIGHT 
39.75E +00 


39.576 +00 


40.09E +00 


soegeonul 


38.35E +00 


38.98E +00 


39.71E+00 


40.09E+00 


38.85E +00 


IN DEGREES Ceseis 


LEFT 
39.90&+00 


38.83E+00 
30.99E*00 
38.096 +00 
38.17E+00 
39.526 400 
39. Soe yuu 
40.20E +00 


39.61E+00 


BOTTOM 
SJ239e200 


S3.o 6 400 


40.59& +00 


39.48E +00 


339.60£+00 


40.48E+00 


39.89&+00 


33.576 °00 


38.95E+00 


eee 
10.04€+00 
10.06& +00 


RIGHT 
10.00€+00 
10.146 +00 


HEAT EXCHANGERS TEMPERATURES: 
BOTTOM: . 


TOP: 


PLANE TEMPERATURES ARE : 


1(SS)>) 25.695 700 
T(56):3 27.76E+00 
-42E +00 
-37E +00 
.87E+00 
-23E +8) 
.04E +00 
~41£ +00 
196 +00°- 


BACK 
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BACK 
45.05E+00 


48 .S8E+00 


46.17E+00 


45.71E+00 


45.28E+00 


S6 .04E+00 


S3E +00 


46. 


Sueesee c00 


49.58E +00 


TEMPERATURE DATA FOR INPUT POWER Qin = 


CHIP NO}: 
POWE 
Gnip N02: 
POWER 


CHIP NQO3: 
POWER 


CHIP NOQ4: 
POWcR 


GmlPe NOS: 
POWE 
CHIP NO6: 
POWER 


CHIP NOQ7: 
POWER 


CHIP NOS: 
POWE 


GHle sNOS ; 
POWE 


AMBIENT TEMP WAS; 
VOLIMETER READING WAS: 
TEMP WAS: 


BATH 


TABLE 7 


THESE RESULTS °ARE SIUREDMI F lie 


conc 


lO. C 


081011130 


5.244 «#61 


THIS RUN WAS RECORDED WITH FLOW VIZ 
REE TEMPERATURE SSoee tl PEGREES CEeLelus 


CENTER 
46.68E+0 


R (WATTS): 


‘Gee oc +00 
(WATTS): 


45.70E+00 
(WATTS): 


46.22E+00 
(WATTS): 


45.4SE +00 


R (WATTS): 


43.SS5E+00 
(WATTS): 


47.36E+00 
(WATTS): 


46.0S5E+00 


R (WATTS): 


44,78E+00 


R (WATTS): 


HEAT EXCHANGERS 
BOTTOM: - 


TOP: 


TOP 


43.58E+00 
27.60E-01 


43.02E+90 
27.63E-01 


43.18&+00 
ereece-Ot 


42.88E+90 
Bi weoe ai 


40.31&+00 
2) 2o7E <0! 


45.08E+00 
e/.29E-01 


44.69E+00 
27.60c-01 


43.926 +00 
27,7 6E= 0) 


ag. eck +00 
Zlso5e 04 


TEMPERATURES: 


RIGHT 
43.37€+00 


43.14E+00 


43.85E+00 


43.55E +00 


41.67E+00 


41.98& +00 


43.43E+00 


43.S56e+00 


41.99&+00 


BACK PLANE TEMPERATURES ARE : 


T(SS): 26 
T(S6): 29 
T(S7):) 31 
3 2a) es 
yay: 2 

TO7ahs 29 
TO75): 30 
T(76): 28 
TS oe 


.87E+00 
.03SE+00 
~ 13E+00 
-93E+00 
ooo +00 
ei teal ie: 
~76E +10 
~82E +00 
.57E +00 
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RIGHT 
ao.c1E-0) 


43 


&2. 


& 


42 


a} 


4] 


43. 


43. 


43. 


42. 


BERT 


69E +00 


04E+00 


» 33E +00 


.566+00 


oor? 00 


23E+00 


20E+00 


77E+00 


90E+00 


LEFT 


BOTTOM 


43.06E +00 
42.866 +00 
44.32E+00 
43.20E +00 
43.18E+00 
44.00E+00 
43.566 +00 
43.346 +00 


41.90E+00 


Soe 99-0) 
10.00£*00 98.46E-01 


2a! 


W 


BACK 


49.816 +00 
S2e2E +00 
50.88E+00 
S0.62E+00 
49.594E+00 
63.01E+00 
S51.8t&+00 
S7.24E+00 


SS.3S5E +00 


TABEE 8 


TEMPERATURE DATA FOR INPUT POWER Q;, = 3.1 W 


THESE RESULTS ARE STORED IN FILE: Da8t011S51S 
AMBIENT TEMP WAS: 23°56 
VOLIMETCR FEADING ee 6.8019 #61 


BATH TEMP AS: 10 
THIS RUN WAS RECORDED WITH FLOW VIZ 


ALL TEMPERATURES ARE IN DEGREES CELCIUS 


CENTER TOP RIGHT LEFT BOTTOM BACK 
CHIP NOI}: SO.27E+00 46.37E6+00 46.67&+00 46.366*00 46.34£+00 S3.20E+00 
POWER (WATTS): 31.S8E-01 
CHIP NO2: 49.47£+00 45.72E6+*00 45.94€+00 44.88E+900 4S5.6S5E+00 S8.34E+00 
POWER (WATTS): 31.60E-01 : 
CHIP NO3: 48.49&€+00 45.7S5E+00 46.70€+00 44.97E+00 47.04E+00 S4.80£+00 
POWER (WATTS): 31.S0E-01 
CHIP NO4: 49.72£*00 46.06€+00 46.89&+00 44.6S5E*00 46.S1E+00 S5.03E+00 
POWER (WATTS): 31.19E-01 
CHIP NOS: 48.678E+00 43.79&+00 44.88E6+00 44.S0E+*00 46.4SE+00 54.10&+00 
POWER (WATTS): 31.32E-01 
CHIP NO6: 46.42E+00 48.26€+00 44.90E+00 46.17E+00 47.08E€+00 68.15E+00 
POWER (WATTS): 31.23E-01 
CHIP NO7: 51.67E+00 48.20€+00 46.82E+00 46.S51E+00 46.95E&+00 56.86E+00 
POWER (WATTS): 31.650€-01 
CHIP NOS: 49.44€+00 46.99&+00 46.61E+*00 46.80E*00 46.32E+00 62.31E+00 
POWER (WATTS): 31.80E-01 
CHIP NO9: 48.20€+00 46.19€+00 44.97E+00 46.00€+00 45.2S5E+00 60.81E+00 
POWER (WATTS): 31.S3E-01 
HEAT EXCHANGERS TEMPERATURES: RIGHT LEFT 
BOTTOM: | 10.02E+00 10.05E+00 
TOP: 10.09E+*00 99.64E-0!1 
BACK PLANE TEMPERATURES ARE : 
T(SS): 28.S8E+00 
T(S6): 30.93E+00 
T($7): 33.10&+00 
T(72): 31.64E+00 
T(73):3 31.018 +00 
1174): 31.47€+00 
T(75): 32.S6E+00 
T(76): 30.46E+00 
T(77)3) =©330.21£ +090 
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TABLE 9 


REDUCED DATA FOR INPUT POWER Qi, = 0.3 W 


THE RAW Emt DATA ARE FROM THE FILE: D80700925 
THE POWER SETTING PER CHIP WAS: 0.3 WATTS 
QNET(H) Tava-Ts Nu ZUNC IN Nu 
27 see -O02 56.03E=01 10.61E£+00 207556 -01 


TEMPERATURE BASED RAYLEIGH NUMBER #* E-6 IS: 104.43E-92 


% UNCERTAINTY IN THE TEMPERATURE BASED RAYLEIGH NUMBER IS :179. 


PEUX BROED RAYLEIGH uNUMBER@@ E-7 IS: > 110.766-02 
Z UNCERTAINTY IN FLUX BASED RAYLEIGH NUMBER IS: 100.88E-02 


Ae oees0e Sle Sal PO 336400 20.14E-01 
TEMPERATURE BASED RAYLEIGH NUMBER * E-6 IS: 107.44E-02 


%Z UNCERTAINTY IN THE TEMPERATURE BASED RAYLEIGH NUMBER IS :174. 


REUXSeGoeD Rey LE LGne NUMBER “9b -/ eto ee 110.98E-02Z 
Z UNCERTAINTY IN FLUX BASED RAYLEIGH NUMBER IS: 100.88E-02 


Ai coe 02 B77 OES I 10.28E+00 20 lige = 01 
TEMPERATURE BASED IRATECIGHINUMBER * E-6 15: 107.63E-02 
4 UNCERTAINTY IN THE TEMPERATURE BASED RAYLEIGH NUMBER IS :174 
neEUX BASEDORAYLELCH NUMBER * E-/7 [Ss "110.63E-02 
Z UNCERTAINTY IN FLUX BASED RAYLEIGH NUMBER IS: 100.88E-02 


GU SUE Ss Ye 55. 156-01 10.65E+00 20.82E-01 
TEMPERATURE BASED RAYLEIGH NUMBER * E-6 IS: 102.58E-02 


% UNCERTAINTY IN THE TEMPERATURE BASED RAYLE1GH NUMBER IS :182. 


Pex BASED RAPE IGHONUMEERS*) E-7 IS; 109.23E-02 
% UNCERTAINTY IN FLUX BASED RAYLEIGH NUMBER IS: 100.89E-02 


2/2 Ge 02 55.S8E-01 10.42E+00 20.42E-01 
TEMPERATURE BASED RAYLEIGH NUMBER * E-6 IS: 105.40E-02 
% UNCERTAINTY IN THE TEMPERATURE BASED RAYLEIGH NUMBER IS :177 
BEUX BASED RAYE IGH NUMBER * &-7.15:  109.36E-02 
% UNCERTAINTY IN FLUX BASED RAYLEIGH NUMBER IS: 100.89E-02 


AT elie 0e Dolor SOT [O213E400 20.00E-01 
TEMPERATURE BASED RAYLEIGH NUMBER * E-6 IS: 108.S54E-02 


Z%Z UNCERTAINTY IN THE TEMPERATURE BASED RAYLEIGH NUMBER IS :172. 


FLUX BASED RAYLEIGH NUMBER * E-7 IS: 109.94E-02 
Z UNCERTAINTY IN FLUX BASED RAYLEIGH NUMBER IS: 100.89E-02 


2h -45b.-02 Sas) Sn 10.62E+00 Z20352E-01 
TEMPERATURE BASED RAYLEIGH NUMBER * E-6 IS: 104.64E-02 


% UNCERTAINTY IN THE TEMPERATURE BASED RAYLEIGH NUMBER IS :178. 


FLUX BASED RAYLEIGH NUMBER * E-7 IS: 111.08E-02 
Z UNCERTAINTY IN FLUX BASED RAYLEIGH NUMBER IS: 100.88E-02 


2BGGE 02: S8.00E-01 10.33£+00 20.04E-01 
TEMPERATURE BASED RAYLEIGH NUMBER * E-6 IS: 108.23E-02 


% UNCERTAINTY IN THE TEMPERATURE BASED RAYLEIGH NUMBER IS :173. 


FEU @epeeo RAYLEIGH NUMBER * E-7 IS: Ti1.7SE-02 
%Z UNCERTAINTY IN FLUX BASED RAYLEIGH NUMBER IS: 100.87E-02 


2) 35e U2 SS eese 0) 10.37E+00 Z20ma236-01 
TEMPERATURE BASED RAYLEIGH NUMBER * E-6 IS: 106.74E-92 


% UNCERTAINTY IN THE TEMPERATURE BASED RAYLEIGH NUMBER IS :175. 


FLUX BASED RAYLEIGH NUMBER # E-7 IS: 110./70E-02 
% UNCERTAINTY IN FLUX BASED RAYLEIGH NUMBER 1S: 100.88E-02 
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Gab=02 


P49E-02 


66E-02 


Ese 


[Sle] 2: 


40E-02 


TABLE 10 


REDUCED DATA FOR INPUT POWER Q:, = 0.7 W 


THE RAW Emf DATA ARE FROM THE FILE: 030701535 
THE POWER SETTING PER CHIP WAS: 0.7 WATTS 
CHIR QNET (CW) Tava-Ts Nu ZUNC IN Nu 
1 (0. sUE {02 PS TES 1Seoce ae 0) lse3se-01 


TEMPERATURE BASED RAYLEIGH NUMBER * E-6 IS: 228.83E-02 

% UNCERTAINTY IN THE TEMPERATURE BASED RAYLEIGH NUMBER IS :872.27E-03 
FLUX BASED RAYLEIGH NUMBER # E-7 IS: 304.72c-02 

Z UNCERTAINTY IN FLUX BASED RAYLEIGH NUMBER IS: 101.0SE-02 


2 70mG2e 02 ~S1E +00 132326400 ee Sor 01 
TEMPERATURE BASED RAYLEIGH NUMBER * E-6& IS: 228.84E-02 
% UNCERTAINTY IN THE TEMPERATURE BASED RAYLEIGH NUMBER IS :872.24E-03 
FLUX BASED RAYLEIGH NUMBER * E-7 IS: 304.80E-02 
%Z UNCERTAINTY IN FLUX BASED RAYLEIGH NUMBER IS: 101.0SE-02 


3 707 09E-02 liso 16200 13.288200 13. 356-01 
TEMPERATURE BASED RAYLEIGH NUMBER * £-6 IS: 228.78E-02 
% UNCERTAINTY IN THE TEMPERATURE BASED RAYLEIGH NUMBER IS :872.44E-03 
FLUX BASED RAYLEIGH NUMBER * E-7 IS: 303.80E-02 
Z UNCERTAINTY IN FLUX BASED RAYLEIGH NUMBER IS: 101.06E-02 


5 69.41E-02 ialecoe +00 13.32E +00 13.43E-01 
TEMPERATURE BASED RAYLEIGH NUMBER * E-6 IS: 225.48E-02 
% UNCERTAINTY IN THE TEMPERATURE BASED RAYLEIGH NUMBER IS :883.80E-03 
FLUX BASED RAYLEIGH NUMBER * E-7 IS: 300.34E-02 
% UNCERTAINTY IN FLUX BASED RAYLEIGH NUMBER IS: 101.07E-02 


S bo. UE a02 ite coet00 13. 936+00 13. 4980 
TEMPERATURE BASED RAYLEIGH NUMBER #* E-6 IS: 222.64E-02 
% UNCERTAINTY IN THE TEMPERATURE BASED RAYLEIGH NUMBER IS :893.83E-03 
FLUX BASED RAYLEIGH NUMBER * E-7 IS: 301.13E-02 
4 UNCERTAINTY IN FLUX BASED RAYLEIGH NUMBER IS: 101.06E-02 


6 621502 11.34E+00 13.38E+00 13.44E-01 
TEMPERATURE BASED RAYLEIGH NUMBER * E-6 IS: 224.97E-02 
%Z UNCERTAINTY IN THE TEMPERATURE BASED RAYLEIGH NUMBER IS :885.S56E-03 
FLUX BASED RAYLEIGH NUMBER * E-7 1S eie2e-02 
7% UNCERTAINTY IN FLUX BASED RAYLEIGH NUMBER IS: 101.07&-02 


7 70.49E-02 VieS7E+00 13. 28E+00 liseoze sO} 
TEMPERATURE -BASED RAYLEIGH NUMBER * E-6 IS: 230.19E-02 
Z UNCERTAINTY IN THE TEMPERATURE BASED RAYLEIGH NUMBER IS :867.69E-03 
FLUX BASED RAYLEIGH NUMBER * E-7 IS: 305./75E-02 
% UNCERTAINTY IN FLUX BASED RAYLEIGH NUMBER IS: 101.0S5E-02 


8 70/7 BER 02 we 1iLS62E+008 13.28E+00 13.30E-01 
TEMPERATURE BASED RAYLEIGH NUMBER * E-6 IS: 231.24E-02 
% UNCERTAINTY IN THE TEMPERATURE BASED RAYLEIGH NUMBER IS :864.18E-03 
FLUX BASED RAYLEIGH NUMBER * E-7 [S: 307.15E-02 
% UNCERTAINTY IN FLUX BASED RAYLEIGH NUMBER IS: 101.05E-02 


9 70m) 7ES02 41.34E+00 13.49E+00 13.43E-01 
TEMPERATURE BASED RAYLEIGH NUMBER * E-6 IS: 225.08E-02 
% UNCERTAINTY IN THE TEMPERATURE BASED RAYLEIGH NUMBER IS :885.19E-03 
FLUX BASED RAYLEIGH NUMBER #* E-7 IS: 303.58E-02 
%Z UNCERTAINTY IN FLUX BASED RAYLEIGH NUMBER IS: 101.06E-02 
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TABLE 11 


REDUCED DATA FOR INPUT POWER Qin = 1.1 W 


THE RAW Ent DATA ARE FROM THE FiILEéE: 080711035 
THE FOWER SETTING PER CHIP HAS: 1.1 WATTS 
QNET (CW) Tava-Ts Nu ZUNC IN Nu 
11.45E-0! 16.676 +00 (S201 +00 fee? 7e 01 


TEMPERATURE BASED RAYLEIGH NUMBER «+ E-6 IS: 343.22£-02 

% UNCERTAINTY IN THE TEMPERATURE BASED RAYLEIGH NUMBER IS :602. 
BEUACROED WAT EelGH NUINBeRe* E-¢ Ig.) 515.31E-O2 

4 UNCERTAINTY IN FLUX BASED RAYLEIGH NUMBER IS: 101.06E-02 


11.46e-01 16.59E +00 len 0 SE +00 niece 01 
TEMPERATURE BASED RAYLEIGH NUMBER * E-S IS: 341.44E-02 
% UNCERTAINTY IN THE TEMPERATURE BASED RAYLEIGH NUMBER IS :605. 
BU ose RAYEE LG NUMBERG E-7 oig280515.066-02 
% UNCERTAINTY IN FLUX BASED RAYLEIGH NUMBER IS: 101.06E-02 


Nl. 42F-01 16.56E+00 H5e06c +00 is JE -01 
TEMPERATURE BASED RAYLEIGH NUMBER * E-6 IS: 340.72E-02 
% UNCERTAINTY IN THE TEMPERATURE BASED RAYLEIGH NUMBER IS :606 
PEUAM DROBO RAYCEIGH NUMBER * E-7 1S: 513.21E-02 
4 UNCERTAINTY IN FLUX BASED RAYLEIGH NUMBER IS: 101.06E-02 


soe soi 16.41E+00 15.05E+00 11.82E-01 
TEMPERATURE BASED RAYLEIGH NUMBER * E-6 IS: 337.05E-02 
4% UNCERTAINTY IN THE TEMPERATURE BASED RAYLEIGH NUMBER IS :612. 
PEUX DASED RA gee [GH NUMBER * E-7 IS; .507.35E-02 
% UNCERTAINTY IN FLUX BASED RAYLEIGH NUMBER IS: 101.07E-02 


Iie SSE =O) 16.05E+00 15.45E+00 lic oe — 01 
JENEER ATURE BASED RAYEEICGH NUMBER * E-6 IS: 328.29E-02 
% UNCERTAINTY IN THE TEMPERATURE BASED RAYLEIGH NUMBER IS :626. 
FLUX BASED RAYLEIGH NUMBER * E-7 IS; 5S07.3S5E-02 
%Z UNCERTAINTY IN FLUX BASED RAYLEIGH NUMBER IS: 101.07E-02 


11.34E-01 lise 22.00 IS so5E +00 ee ogZze 01 
(ehiceh PUR eee AoEDORATEEIGH NUMBER © E-6 153) 9325.25E-02 
% UNCERTAINTY IN THe TEMPERATURE BASED RAYLEIGH NUMBER IS :631. 
FEU SSAoeDeRATEEIGH NUMBER * E=7 IS: 505.88E-02 
% UNCERTAINTY IN FLUX BASED RAYLEIGH NUMBER IS: 101.07E-02 


le 4SE-01 16.81E+00 14.93E+00 lie 7 GE=01 
TEMPERATURE BASED RAYLEIGH NUMBER * E-6 IS: 346.68E-02 
% UNCERTAINTY IN THE TEMPERATURE BASED RAYLEIGH NUMBER IS :597. 
FLUX BASED RAYLEIGH NUMBER * E-7 IS: 517.48E-02 
% UNCERTAINTY IN FLUX BASED RAYLEIGH NUMBER IS: 191.06E-02 


WEES3E-01 16.76E+00 15.03E+00 le SE = 01 
TEMPERATURE BASED RAYLEIGH NUMBER * E-6 IS: 345.42E-02 
% UNCERTAINTY IN THE TEMPERATURE BASED RAYLEIGH NUMBER IS :599. 
BEURee Hoe Oe RAYEEICH NUMBER SE -7 [52 05ol9. 33E-02 
% UNCERTAINTY IN FLUX BASED RAYLEIGH NUMBER IS: 101.05E-02 


ie 436-01 ~- 16.28E+00 15.34E+00 11.84E-01 
TEMPERATURE BASED RAYLEIGH NUMBER * E-6 IS: 333.97E-02 
% UNCERTAINTY IN THE TEMPERATURE BASED RAYLEIGH NUMBER IS :617. 
FLUX BASED RAYLEIGH NUMBER # E-7 IS; 512.18E-02 
%Z UNCERTAINTY IN FLUX BASED RAYLEIGH NUMBER IS: 101.06E-02 
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TABLE 12 


REDUCED DATA FOR INPUT POWER Qj, = 1.5 W 


THE RAW Emf BATA ARE FROM THE FILE: D80720730 
THE POWER SETTING PER CHIP WAS: 1.5 WATTS 
QNET(W) Tava-Ts Nu ZUNC IN Nu 
14. 596-0] 19. 056700 16.76E+00 11.41E-01 


TEMPERATURE BASED RAYLEIGH NUMBER * £-6 [5: 401.43E-02 

Z UNCERTAINTY IN THE TEMPERATURE BASED RAYLEIGH NUMBER IS :527 
FLUX BASED RAYLEIGH NUMBER * E=7 1S Bib7 2s) B02 

% UNCERTAINTY IN FLUX BASED RAYLEIGH NUMBER IS: 10%.12E-02 


14.59E-01 19.08E+00 16.74E+00 11.40E-01 
TEMPERATURE BASED RAYLEIGH NUMBER * E-6 IS: 402.2SE-02 


7% UNCERTAINTY IN THE TEMPERATURE BASED RAYLEIGH NUMBER IS :526. 


FLUX BASED RAYLEIGH NUMBER * Es? IS .eG7 sai 02 
% UNCERTAINTY IN FLUX BASED RAYLEIGH NUMBER IS: 101.12E-02 


14.54E-01 iS3gSe<00 16.45E+00 Mie ove U1 
TEMPERATURE BASED RAYLEIGH NUMBER * E-6 IS: 409.23E-02 
% UNCERTAINTY IN THE TEMPERATURE BASED RAYLEIGH NUMBER IS :519 
FLUX BASED RAYLEIGH NUMBER * E-7 IS: 673.03E-02 
%Z UNCERTAINTY IN FLUX BASED RAYLEIGH NUMBER IS: 101.12E-02 


-40E-01 18. 73E+00 16.83E+00 11.45E-01 
TEMPERATURE BASED RAYLEIGH NUMBER * B-6 (S-9 39a .a5E-02 


7% UNCERTAINTY IN THE TEMPERATURE BASED RAYLEIGH NUMBER IS :536. 


FLUX BASED RAYLEIGH NUMBER * E=7 IS7eeS61 .82E-02 
% UNCERTAINTY IN FLUX BASED RAYLEIGH NUMBER IS: 101.13€-02 


14.46E-01 18.39E+00 iiecue 00 is49e-01 
TEMPERATURE BASED RAYLEIGH NUMBER * E-6 IS: 385.00E-02 


4 UNCERTAINTY IN THE TEMPERATURE BASED RAYLEIGH NUMBER I[S :546. 


FLUX BASED RAYLEIGH NUMBERS 6-7 915. beze05e -02 
%Z UNCERTAINTY IN FLUX BASED RAYLEIGH NUMBER IS: 101.13&-02 


14.44E-01 18.59E +00 16.99£+00 11.47E-01 
TEMPERATURE BASED RAYLEIGH NUMBER * E-6 IS: 389.86E-02 


% UNCERTAINTY IN THE TEMPERATURE BASED RAYLEIGH NUMBER IS :540. 


FLUX BASED RAYLEIGH NUMBER * E-7 IS: 662.46E-02 
% UNCERTAINTY IN FLUX BASED RAYLEIGH NUMBER IS: 101.13E-02 


14.62E-01 Wee 200 16.60E+00 li coe 0 | 
TEMPERATURE BASED RAYLEIGH NUMBER #* E-6 IS: 407.23E-02 
% UNCERTAINTY IN THE TEMPERATURE BASED RAYLEIGH NUMBER IS :521 
FLUX BASED RAYLEIGH NUMBER * E-7 [Seep Ue 
%Z UNCERTAINTY IN FLUX BASED RAYLEIGH NUMBER IS: 101.11E-02 


14.68E£-01 — 19, 486 +00 16.S0E+00 11.3S5E-01 
TEMPERATURE BASED RAYLEIGH NUMBER * E-6 IS: 


A 2 5306202 
% UNCERTAINTY IN THE TEMPERATURE BASED RAYLEIGH NUMBER IS :516. 


FLUX BASED RAYLEIGH NUMBER * E-7 IS: 680.47E-02 
% UNCERTAINTY IN FLUX BASED RAYLEIGH NUMBER IS: 101.11E-02 


14.S56E-01 *'19.09E+00 16.69£+00 11.40E-01 
TEMPERATURE BASED RAYLEIGH NUMBER * E-6 IS: 402.60E-02 


% UNCERTAINTY IN THE TEMPERATURE BASED RAYLEIGH NUMBER IS :526. 


FLUX BASED RAYLEIGH NUMBER * E-7 IS: 671.79E-02 
% UNCERTAINTY IN FLUX BASED RAYLEIGH NUMBER IS: 101.12E-02 
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REDUCED DATA FOR INPUT POWER Orn = 1.9 W 


THE RAW Emt DATA ARE FROM THE FILE: 089721819 
THE POWER SETTING PER CHIP WAS: led Halts 
QNET (W) Tava-is Nu ZUNC IN Nu 
io aage-0 1 eo ce t0) aoe 010 10.88 -01 


TEMPERATURE BASED RAYLEIGH NUMBER + E- PE tS oo on Oe 

% UNCERTAINTY IN THE TEMPOCRATURE BASED RAYLEIGH NUMBER IS :403 
FLUX BASED RAYLEIGH NUMBER « E-7 IS: 947.69E-02 

4 UNCERTAINTY IN FLUX BASED RAYLEIGH NUMBER IS: 101.07E-02 


13 00E=0;] 24.72E+00 62832 +00 LOee je = OT 
TEMPERATURE BASED RAYLEIGH NUMBER * E-6 IS: S60.03E-02 


“% UNCERTAINTY IN THE TEMPERATURE BASED RAYLEIGH NUMBER [S :406. 


MEU ARSED RATEEIGH NUMBERe* E-7 JS:) 945.98E-C2 
%Z UNCERTAINTY IN FLUX BASED RAYLEIGH NUMBER IS: 101.07E-02 


18.94E-01 24.68E+00 16.866 +00 He o0E-0% 
TEMPERATURE BASED RAYLEIGH NUMBER * E-6 IS: 5S59.04E-02 


%Z UNCERTAINTY IN THE TEMPERATURE BASED RAYLEIGH NUMBER IS :40/7. 


PEUX Baseoenarleied NUMBER * &-/7 IS:> 942.S2E-02 
4 UNCERTAINTY IN FLUX BASED RAYLEIGH NUMBER IS: 101.08E-02 


e275 -01 24.53E +00 16.80E+00 Omni = 0 } 
NEGEeERAIURE BASED RAYEEIGH NUMBER * E=6 IS: 554.73E-02 
% UNCERTAINTY IN THE TEMPERATURE BASED RAYLEIGH NUMBER IS :409 
PUUX ShoeOeRaVeelGHe NUMBER #7 E=7 IS: 931./79E-02 
%Z UNCERTAINTY IN FLUX BASED RAYLEIGH NUMBER IS: 101.09E-02 


ieeooc 70! 23.88E+00 iy nocet00 Oe Soe 0 } 
TEMPERATURE BASED RAYLEIGH NUMBER # E-6 IS: 5$36.40E-02 
% UNCERTAINTY IN THE TEMPERATURE BASED RAYLEIGH NUMBER IS :420 
PEOASEAoeD RonNLeIGH NUMBER * E-7 152) 928.99E-02 
Z UNCERTAINTY IN FLUX BASED RAYLEIGH NUMBER IS: 101.08E-02 


18.80E-01 c3.83E+00 lieeoces00 ie 25e=0 | 
TEMPERATURE BASED RAYLEIGH NUMBER * E-6 IS: 535.00E-02 
Z% UNCERTAINTY IN THE TEMPERATURE BASED RAYLEIGH NUMBER IS :421 
FLUX BASED RAYLEIGH NUMBER * E-7 IS: 926.84E-02 
%Z UNCERTAINTY IN FLUX BASED RAYLEIGH NUMBER IS: 101.09E-02 


19.04E-01 25.23E+00 Homo Je +00 10.86E-01 
TEMPERATURE: BASED RAYLEIGH NUMBER * E-6 IS: 574.77E-02 
% UNCERTAINTY IN THE TEMPERATURE BASED RAYLEIGH NUMBER IS :398 
GeUpelseD RAYLEIGH NUSBER * E=/ -[S: 953.30E-02 
%Z UNCERTAINTY IN FLUX BASED RAYLEIGH NUMBER IS: 101.07E-02 


S26 =O). 25.29E+00 16.68E+00 10866 -01 
TEMPERATURE BASED RAYLEIGH NUMBER * E-6 IS: 573.88E-02 


% UNCERTAINTY IN THE TEMPERATURE BASED RAYLEIGH NUMBER IS :398. 


FLUX BASED RAYLEIGH NUMBER + E-7 IS: 957.00E-02 
% UNCERTAINTY IN FLUX BASED RAYLEIGH NUMBER IS: 101.07E-02 


18.96E-01 -24.47E+00 17.02E+00 HOmahe- Ot 
TEMPERATURE BASED RAYLEIGH NUMBER * E-6 IS: 552.91E-02 


% UNCERTAINTY IN THE TEMPERATURE BASED RAYLEIGH NUMBER IS :410. 


FLUX BASED RAYLEIGH NUMBER * E-7 IS: 941.24E-02 
% UNCERTAINTY IN FLUX BASED RAYLEIGH NUMBER IS: 101.08E-02 


5S 


wie=03 
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TABLE 14 


REDUCED DATA FOR INPUT POWER Qi, = 2.3 W 


THE RAW Emr DATA ARE FROM THE FILE: D80731030 
THE POWER SETTING PER CHIP WAS: 2.3 WATTS 
CHIP QNET(W) Tavg-Ts Nu ZUNC IN Nu 
1 2c.8S5E-01 23.31E+00 17.76E+00 10.71E-01 


TEMPERATURE BASED RAYLEIGH NUMBER * E-& IS: 658.S53E-02 

% UNCERTAINTY IN THE TEMPERATURE BASED RAYLEIGH NUMBER IS :355.02E-03 
FLUX BASED RAYLEIGH NUMBER * £-7 IS: 116.98E-01 

Z UNCERTAINTY IN FLUX BASED RAYLEIGH NUMBER IS: 101.07E-02 


2 22.86E-01 Cie OU 17.98E+00 ieee =0t 
TEMPERATURE BASED RAYLEIGH NUMBER * E-6 IS: 648.40E-02 
% UNCERTAINTY IN THE TEMPERATURE BASED RAYLEIGH NUMBER IS :359.32E-03 
FLUX BASED RAYLEIGH NUMBER + E-7 IS: 116.60E-01 
Z%Z UNCERTAINTY IN FLUX BASED RAYLEIGH NUMBER IS: 101.07E-02 


3 Zee GeO af oe +00 mf GO +00 non 7 SE-01 
TEMPERATURE BASED RAYLEIGH NUMBER + E-6 IS: 648.36E-02 
~Z UNCERTAINTY IN THE TEMPERATURE BASED RAYLEIGH NUMBER IS :359.34E-03 
FLUX BASED RAYLEIGH NUMBER 2eE—7 01S. Seigeecceaon 
% UNCERTAINTY IN FLUX BASED RAYLEIGH NUMBER IS: 101.07E-02 


4 e2esoe 01 2¢/.37E+00 17 so 1E+00 10.73E-01 
TEMPERATURE BASED RAYLEIGH NUMBER * E-6 IS: 645.46E-02 
% UNCERTAINTY IN THE TEMPERATURE BASED RAYLEIGH NUMBER IS :360.59E-03 
FLUX BASED RAYLEIGH NUMBER * E-7 IS: 114.97E-01 
Z UNCERTAINTY IN FLUX BASED RAYLEIGH NUMBER IS: 101.08E-02 


5 Censor su) 27 oes e 18.38E+00 10.77E-01 
TEMPERATURE BASED RAYLEIGH NUMBER *« E-6 IS: 622.74E-02 
Z UNCERTAINTY IN THE TEMPERATURE BASED RAYLEIGH NUMBER IS :370.80E-03 
FLUX BASED RAYLEIGH NUMBER * E-7 IS: 114.48E-01 
Z UNCERTAINTY IN FLUX BASED RAYLEIGH NUMBER IS: 101.08E-02 


6 €2-6UE=01 cfs (SE+00 18.31E+00 10.76E-01 
TEMPERATURE BASED RAYLEIGH NUMBER * £-6 IS: 624.26E-02 
%, UNCERTAINTY IN THE TEMPERATURE BASED RAYLEIGH NUMBER IS :370.10E-03 
FLUX BASED RAYLEIGH NUMBER * E-7 [S: 114.28E-01 
% UNCERTAINTY IN FLUX BASED RAYLEIGH NUMBER IS: 101.08E-02 


7 22.90E-01 28.72E+00 17 .55E+00 10.70E-01 
TEMPERATURE BASED RAYLEIGH NUMBER * E-6 IS: 671.04E-02 
Z UNCERTAINTY IN THE TEMPERATURE BASED RAYLEIGH NUMBER IS :349.89E-03 
FLUX BASED RAYLEIGH NUMBER * E-7 IS: 117.80E-01 
% UNCERTAINTY IN FLUX BASED RAYLEIGH NUMBER IS: 101.07E-02 


8 23.00E-01 — 20 soge 500 17.71E+00 10.70E-01 
TEMPERATURE BASED RAYLEIGH NUMBER * E-6 IS: 667.0SE-02 
% UNCERTAINTY IN THE TEMPERATURE BASED RAYLEIGH NUMBER IS :351.51E-03 
FLUX BASED RAYLEIGH NUMBER * E-7 IS: 118.13E-01 
% UNCERTAINTY IN FLUX BASED RAYLEIGH NUMBER IS: 101.06E-02 


eS) 22 ec0E-0)) *' 27.68E+00 18.1sE+00 10.74E-01 
TEMPERATURE BASED RAYLEIGH NUMBER * E-6 IS: 639.69E-02 
% UNCERTAINTY IN THE TEMPERATURE BASED RAYLEIGH NUMBER IS :363.12E-03 
FLUX BASED RAYLEIGH NUMBER #* E-/ IS: 115.97E-01 
%Z UNCERTAINTY IN FLUX BASED RAYLEIGH NUMBER IS: 101.07E-02 
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TAS bE S 


REDUCED DATA FOR INPUT POWER Qj, = 2.7 W 


TRE RAW Emt DATA ARE FROM THE “ILE: Ne0731755 
THE PCHER SETTING PER CHIP WAS: 2.7? wATTS 
CHIP QNET(W) Tava-Ts Nu ZUNC IN Nu 
1 PETER GIT A Ses Oa 33.28E +090 18.16£+00 ie eo5e— 01) 


fevmesnilune BASED RAYLEIGH NUMBER «5—-6 ISt = 613. 30B-U2 

% UNCERTAINTY IN THE TEMPERATURE BASED RAYLEIGH NUMBER IS :302.08c-03 
Plu BASED RAYLEIGH NUMBER * &-7 1S: 14/7./0E-01 

Z UNCERTAINTY IN FLUX .BASED RAYLEIGH NUMBER IS: 101.05E-02 


2 Mix SIE = 01 Se? 7 E+00 18.44E +00 EORSoE- 01 
TEMPERATURE BASED RAYLEIGH NUMBER * E-6 I5: 796.9SE-02 
% UNCERTAINTY IN THE TEMPERATURE BASED RAYLEIGH NUMBER IS :306.72E-03 
FLUX BASED RAYLEIGH NUMBER # E-7 IS: 146.98E-01 
% UNCERTAINTY IN FLUX BASED RAYLEIGH NUMBER IS: 101.0S5E-02 


3 27e30E -01 22 .55E4U0 187 o25-00 HORS 7E-01 
TEMPERATURE BASED TRAYEETGH NUMBER * E=-6 [5: 7389.28E-02 
% UNCERTAINTY IN THE TEMPERATURE BASED RAYLEIGH NUMBER IS :308.97E-03 
PEUGERoeD RATCEICH NUMBER * E-7 15: 146.13E-07 
4 UNCERTAINTY IN FLUX BASED RAYLEIGH NUMBER IS: 101.0SE-02 


4 21.0320 92.756 +00 16.226 +400 Or OE= 0h 
TEMPERATURE BASED RAYLEIGH NUMBER * E-6 IS: 796.28E-02 
% UNCERTAINTY IN THE TEMPERATURE BASED RAYLEIGH NUMBER IS :306.91E-03 
FLUX BASED RAYLEIGH NUMBER * E-7 IS: 145.04E-01 
% UNCERTAINTY IN FLUX BASED RAYLEIGH NUMBER IS: 101.06E-02 


S ZT eVSe- Ot ay. 7 3E 700 18.87E+00 10.S53E-01 
TEMPERATURE BASED RAYLEIGH NUMBER + E-6 IS: 763.S1E-02 
%Z UNCERTAINTY IN THE TEMPERATURE BASED RAYLEIGH NUMBER IS :316.79E-03 
FLUX BASED RAYLEIGH NUMBER #* E-7 IS: 144.06E-01 
% UNCERTAINTY IN FLUX BASED RAYLEIGH NUMBER IS: 101.06E-02 


6 2707-0) ol 7de-00 16.78E+00 FOSSIE-01 
TEMPERATURE BASED RAYLEIGH NUMBER * E-6 IS: 765.48E-02 
% UNCERTAINTY IN THE TEMPERATURE BASED RAYLEIGH NUMBER IS :316.18E-03 
FLUX BASED RAYLEIGH NUMBER * E-7 IS: 143.75€-01 
% UNCERTAINTY IN FLUX BASED RAYLEIGH NUMBER IS: 101.06E-02 


7 27.44E-01 336e2e+00 17.92E +00 10. S3E—0" 
TeneeRAlLURE BASED RAYLEIGH NUMBER * E-5 [S:> 837. 14E-02 
% UNCERTAINTY IN THE TEMPERATURE BASED RAYLEIGH NUMBER IS :297.22E-03 
FLUX BASED RAYLEIGH NUMBER * E-7 IS: 148.91E-01 
Z UNCERTAINTY IN FLUX BASED RAYLEIGH NUMBER IS: 101.0SE-02 


8 2/2566-01 33.41E+00 VSa2Te +00 10.S4E-01 
TEMPERATURE BASED RAYLEIGH NUMBER * E-6 IS: 817.S56E-02 
% UNCERTAINTY IN THE TEMPERATURE BASED RAYLEIGH NUMBER IS :300.90E-03 
FLUX BASED RAYLEIGH NUMBER * E-7 IS: 148.90E-01 
% UNCERTAINTY IN FLUX BASED RAYLEIGH NUMBER IS: 101.04E-02 


5 Pb fe SAE), OZ 71 IE +00 18.77E+00 TO oSe-01 
TEMPERATURE BASED RAYLEIGH NUMBER * E-6 IS: 775.54E-02 
% UNCERTAINTY IN THE TEMPERATURE BASED RAYLEIGH NUMBER IS :313.07€-03 
FLUX BASED RAYLEIGH NUMBER * E-7 IS: 145.S9E-01 
%Z UNCERTAINTY IN FLUX BASED RAYLEIGH NUMBER IS: 101.0S5E-02 


SW 


1. Graphical Representation of Reduced Data 


The experimental surface temperature data were 
plotted as a dimensionless Nusselt Number vs. both a 
temperature based and a flux based Rayleigh Number, Figures 
4.15 through 4.32. 

Representative error bounds were plotted in Figures 
4.15 and 4.16 for input powers of 1.5 and 1.9 watts. All 
points, with the exception of 1.5 watts were well within the 
error bounds. The point representing 1.5 watts was either 
too high on the Nusselt Number axis or too low on the 
Rayleigh number axis. 

The factors affecting the Nusselt Number were h, L & 
K. Of these three, only h was calculated with uncertainty 
in this experiment. The factors affecting the Rayleigh Nun- 
bers were a, L, AT, Qnet, 8, vy and K. Of these only AT and 
Qnet had uncertainties. 

Examination of Tables 4 and 12 resulted in no data 
which appeared to be extreme in either the raw or reduced 
data sets. One possible explanation for this point is the 
changing flow patterns around this power input. This may be 
a transition region from somewhat orderly flow to completely 
chaotic flow. It is recommended that further investigation 
of this region of power settings be conducted. 

The equations for the best fit lines are very close 
to each other suggesting that an overall correlation that 


would produce good approximations is possible. 
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Figure 4.15 Temperature Based Rayleigh Numbers, Chip 1 
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Figure 4.16 Flux Based Rayleigh Numbers, Chip 1 
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Figure 4.17 Temperature Based Rayleigh Numbers, Chip 2 
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Figure 4.18 Flux Based Rayleigh Numbers, Chip 2 
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Figure 4.19 Temperature Based Rayleigh Numbers, Chip 3 
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Figure 4.20 Flux Based Rayleigh Numbers, Chip 3 
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Figure 4.21 Temperature Based Rayleigh Numbers, Chip 4 
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Figure 4.22 Flux Based Rayleigh Numbers, Chip 4 
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Figure 4.23 Temperature Based Rayleigh Numbers, Chip 5 
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Figure 4.24 Flux Based Rayleigh Numbers, Chip 5 
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Figure 4.25 Temperature Based Rayleigh Numbers, Chip 6 
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Figure 4.26 Flux Based Rayleigh Numbers, Chip 6 
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Figure 4.27 Temperature Based Rayleigh Numbers, Chip 7 
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Figure 4.28 Flux Based Rayleigh Numbers, Chip 7 
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Figure 4.29 Temperature Based Rayleigh Numbers, Chip 8 
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Figure 4.30 Flux Based Rayleigh Numbers, Chip 8 
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Figure 4.31 Temperature Based Rayleigh Numbers, Chip 9 
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Figure 4.32 Flux Based Rayleigh Numbers, Chip 9 
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Figure 4.33 
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Figure 4.34 
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Figure 4.35 Temperature Variation for Input Power 1.1 W 
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yah 


(928e@) SwTt, 


oo02 ost ost owt o2t oot os og Ov 
 aaale 


(988) Sowt, 


oo02 oOSBt Ost Ort O2t oot 08 og Ov 
ee ee al ele SS 


ES ES ee) es a ee ny meee eee 


O02 


Oe 


m9 “OE 


@o 


"SE 


QE 


2 °OE 


v “OE 


(23> Ssunzouedwe, 


8 QE 


ZE 


SE 


12 “OE 


= 
© 
m 


@ 0 
oO oO 
m m 
(D> SsnrzouBdueay 


Ze 


Sa) Ob 4 Coa 


L____3]) [ollos) Cog 


Vile tn 


(208) owt), 


002 oset 


ost oOvt oOet oot 
ee eres emer 1 | 


els) 


oe 


Ov 


Oe 


~> 


oo 


"VE 


z N 
i" in 
m m 
(D> Sunaocveduey, 


O 
in 
m 


Temperature Variation for Input Power 1.9 W 


Figure 4.37 
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Temperature Variation for Input Power 2.7 W 
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Figure 4.40 Temperature Variation for Input Power 3.1 W 
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C. TEMPERATURE VARIATIONS 

The temperature variations with time were investigated 
and reported in Figures 4.33 through 4.40. These results 
were motivated by the studies conducted by Liu, Yang and 
Kelleher {Liu et al., 1987, pp. 486-97]. Three thermocouple 
channels, numbers 1, 13 and 32 were scanned individually at 
a rate of 3 times per second for a period of 200 seconds. 
These quantitative results were then plotted. The results 
clearly demonstrate a multiple frequency content signal for 
most of the power levels investigated in this experiment. 

An observation of this experiment was that the output of 
thermocouple #1 remained steady with time throughout the 
increase in power levels. The temperature level for this 
thermocouple was also generally higher than for the other 
two thermocouples sampled. 

A possible explanation for these results is that the 
cold, stagnant fluid layer created by the bottom heat 
exchanger prevented the penetration of warmer fluid below 
the bottom row of components. This is also supported by 
flow visualization. 

One exception is at 3.1 watts. As seen in Figure 4.40, 
there is a large fluctuation of temperature with time. A 
possible explanation for this fluctuation may be that the 
increased energy acquired by the fluid at this power setting 
enabled the flow to force itself into the stagnant region at 


the bottom of the chamber. 
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V. RECOMMENDATIONS 


There 1S a great deal more experimentation. to be 
accomplished with this component configuration. Three 
possible variations to the chamber are : 


- Placing a spacer in the chamber to reduce the chamber 
width. This would act in the same manner as a shrouding 
wall and would affect the heat transfer characteristics. 


- Setting the bottom heat exchanger to a temperature other 
than 10°C. The energy expended to maintain the bottom 
heat exchanger at 10°C is considerable. If reasonable 
components temperatures can be realized by maintaining 
the bottom heat exchanger at a higher temperature level 
then the energy consumption of the system would be 
reduced. 


- Powering different configurations other than the entire 
array. This could possibly offer some explanation of the 
flow patterns observed in the chamber. 

Other possible improvements to software and hardware 
are: 

- Producing a plotting program that would read the data 
Girectly from the data acquisition unit for plots of 
Nusselt vs. Rayleigh number. 

- Implementation of measurement devices that have greater 
sensitivity to minor fluctuations in temperature 
variation. 


- Development of a program that will perform frequency 
analysis on the data recorded by the program Fastscan. 
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APPENDIX A 


SAMPLE CALCULATIONS 


Throughout this appendix, the sample calculations will 
be based on component number one with a power input of 0.71 
watts. 


1. Determination of Input Power 





Figure A.1 Electrical Network of Power Input 


From Equation (3.4): 


Qin = ((3.218 - 0.5278) (0.5278) )/2 
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2. Estimation of Conduction Loss 

Equation 3.3 requires Tg, Tp and Roop. Te and Tp were 
read directly from thermocouples placed immediately behind 
the heaters and on the back of the board, respectively. 
Reon Was calculated by using the constriction resistance 
formulas developed by Mikic and repeated in [Kraus and Bar- 
Cohen, 1983, pp. 76-9]. The conduction losses for components 
3, 5 and 7 were then averaged to yield a total conduction 
loss. 
3. Calculation of Nusselt and Rayleigh Numbers 

The calculations of the convection coefficient from 


Equation (3.5) was: 


Qnet = HA(Tavg ~ Tc) 


h Qnet 
A(T og 7 Ty) 


703 
els 2 al) =) (OLIN) 


= 106.04 W/m? OK 


From Equation (3.10): 


Ay + T 
_ av Cc 
Teiim = 


2 lee Det ed) 
7 2 


Tus, 


Team =e e7 


Using the relationships of [3-M Manual, 1985, pp. 9-60], the 


physical properties of FC-75 were: 


ke = 


O 
l 


< 
ll 


< 
ll 


Les) 
II 


(.65 - 7.8947 x°alo” + Nae aco) 


0.0637 W/m°c 


1000(1.825 - 0.00246 * Teijn) 


1823.21 Kg/m? 


4180(.241111 + 3.7037 x 105% * T>-j,) 
1032.23 J/Kg°C 
107©(1.4074 - 2.964 x 1072 * Tria 

+ 3.8018 x 1074 * Tr), 


= 2.7308 x 100° Sites 


+ 8.1679 x 1079 * Tein) 


1.024 x 107© m2/sec 


.00246 


- * 
NE FPA) .00246 Tesim 


le 37 To eon 
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iaen from Equation (3.1): 


Nu 


Pr 


Grt 


106.04(.008) 
.0637 


DiS Z 


Oc 


10) 


e063 7 
emer 2 P1032 2 


3.38 x 1078 m/sec 


au 
O 


LGA. se ales 


aie se Tee 


ASEM = lee.) COGS Cae 


(1.024 x 107°)? 


7.592 x 104 
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From Equation (3.11), the 


Was: 


Rat = Gry * Pr 


2S ex 106 


temperature based Rayleigh number 


From Equation (3.12), the flux based Rayleigh number was: 


Rag = 3.06 x 10/ 
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APPENDIX B 


UNCERTAINTY CALCULATIONS 


Uncertainty calculations were evaluated using the root 
mean square formulation described in [Beckwith, Buck and 
Marangoni, 1982, pp. 261-292]. The following calculations 


are for components number one at 0.7 watts (see Table 4.10): 


TABLE 16 


UNCERTAINTIES IN VARIABLES 


Variables Uncertainties Reference 
L 107° m Resolution 
dh 0.1 % Thermocouple readings 
Rp 0.02 1% Variation 
Kp hg ees (Touloukian, 1986, p. 970] 
V O00 1Ey Resolution 


The uncertainty in the net power added to the fluid was: 


2 


eZ 
z SQ, oss 


$0. = (SPower ) 


net 


where 
6 Power = 0.01 * Power 
= .0071 W 


&3 


SQloss = 90-955 * Qioss 


= .000159 W 
Therefore, 


Sores SOO gT 


The uncertainty for the other variables used were: 


2 2 
STeiim = (6Tsink + ae 
= (.01% + 0.12)1/2 
= .1414 K 


2 
6Tpeit = OTavg + 9 Sune) oo 
= (0.12 + .012)1/2 


= .1414 K 
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6 Ra 
Ra 
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AEPENDEX “¢ 


SOFTWARE 


Acquire 


! FILE ACQUIRE 
| EDITED BY LT T J BENEDICT, USN 
! READ FILE *READ_ME" 


COM /Co/ D(7) 
DIA Emf (76) ,Power(9),1(76) 


ICORRELATION FACTORS TO CONVERT EMF TO DEGREES CELCIUS 
DATA 0.10086091 .25727.9. -767345.8. 78025596, 
DATA -9247486589,6.98E11,-2.66£13,3.94E1'4 


READ D(+#) 
Rp=2.000 
PRINTER IS 701 
BEEP 

j 


INPUT "ENTER THE INPUT MODE: O*SYS, «FILE, Im 


IF Ime! THEN 
BEEP 
LHPUT "ENTER THE NAME OF THE FILE TO BE READ” .Oldfiles 


PRINT USING "ISX, "" THESE RESULTS ARE STORED IN FILE : °",10A"°:O!dfiles 


ELSE 

BEEP 

INPUT "ENTER THE NAME OF THE NEW FILE” .Newf iles 

PRINT USING °10X,°*° THESE RESULTS ARE STORED IN FILE: °° ,10A°:Newfile$ 
END IF 

PRINT 

PRINT 

' 


INPUT "ENTER THE AMBIANT TEMP’, Temps 

PRINT USING °15X,°"AMBIENT TEMP HAS: °°, 10A°:Temp$ 
INPUT "FLOW VIZ? Y/N°,Ans$ 

INPUT "VOLTMETER READING’.VS 

PRINT USING *15X,"*VOLTMETER READING HAS: °°, 10A";VS 


INPUT "ENTER THE BATH TEMP” ,.BS 
PRINT USING °15X,°"BATH TEMP WAS 7"",10A°°BS 
' 


IF Ans$="Y" THEN PRINT USING "1SX,°°THIS RUN WAS RECORDED WITH FLOW VIZ"", 


IF Im=1 THEN ASSIGN 8File TO Oldfiles 
PRINT 

} 

IF Im=Q THEN 

CREATE BDAT Newfiles,S 

ASSIGN @File 10 Newfiles 

END IF 

t 


| READ DATA 


j 
IF Im*O0 THEN 
OUTPUT 709:°AR AFOO AL79° 
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PRINT USING "14X,°°ALL TEMPERATURES ARE IN DEGREES CELCIUS*** 


TOP RIGHT LEFT BOTTON 


S40 QCUIPUT 7223°F1 Ri Tioeeco: 

S50 ' 

560 FOR f*0 10 76 

570 OUTPUT 709:°AS° 

580 WAIT 1 

S930 ENTER 722:Emf(T) 

600 BEEP 

610 NEXT J 

620 | 

639 OUTPUT ®File:Emf(e) 

640 ! 

650 ELSE 

660 ENTER @File:Emf(«) 

670 END IF 

680 ! 

690 OUTPUT 709;°TD® 

709 ' 

710 FOR I=*0 TO 60 

720 Sum=U 

730 FOR Jel TO 7 

740 Sum*Sum+D( J) emf (1) ° J 

750 NEXT J 

760 T(1)*Sum 

770 HEXT I 

780 ' 

790 FOR I-71 TO 76 

800 Sum*Q 

810 FOR Je! TO 7 

820 Sum=Sum+D( J) #Emi (1) ° J 

830 NEXT J 

840 TOL) Sum 

050 NEXT I 

860 ! 

870 ! POWER CALCULATIONS 

880 ! 

890 Je! 

900 Volt-Emf (61) 

910 ' 

920 FOR I-62 TO 70 

930 Power(J)-Emf (I) #(Volt-Emf (I) )/Rp 

940 JeJ+l 

350 NEXT I 

960 1 

970 BEEP 

980 BEEF 

990 

1000 '! 

1010 PRINT 

1020 =36! 

1030 PRINT USING "12X,°°CENTER 
BACK"*° 

1040 PRINT 

1050 


060 
61 


80 


1081 


0 
0 
TOTO). TON 
0 
0 


Lt ah USING "1X,°°CHIP NOT: °°,S(Z7Z.0DE.2X).2K.Z2Z.DDE° 210). 701), Tle), FC 
0 104), 


PRINT USING °SX,"°POWER (WATTS): °°,ZZ.0DE°:Power(1) 


PRINT 


70 PRINT USING “1tX,"°CHIP NO2: °°,S5(72.DDE.2X).2X,722.DDE°. F(6), 707), 718), 7 


) 
FRINT USING °SX,°°POWER (WATTS): °°,7Z7.DDE°:Power(2) 


PRINT 
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0 


PRINT USING 


Peart), 1117) 


PRINT USING 
PRINT 


PRINT USING 


om eae, | N23} 


PRINT USING 
PRINT 
PRINT USING 


aia Ti 29) 


PRINT USING 
PRINT 
PRINT USING 


~1(34),7(35) 


PRINT USING 
PRINT 
PRINT USING 


~1(40),7T(41) 


PRINT USING 
PRINT 
PRINT USING 


~1(46).7(47) 


PRINT USING 
PRINT 
PRINT USING 


ml tS2 0.52) 


oe USING 


PRINT 
PRINT 


PRINT USING 
PRINT USING 
PRINT USING 
PRINT 

PRINT 

{ 


PRINT USING 
PRINT 

PRINT USING 
PRINT USING 
FRINT USING 
PRINT USIN 
PRINT USING 
PRINT USING 
PRINT USING 
PRINT USING 
PRINT USING 
BEER 
PRINTER IS 1 
i 


Coelxenme te 72hts sek. 022.00E.2X)"s 


Tk oceeriem Uc mmensic DDE vex) maxene, O0Es leit 13).1(14) 
"SX."°POWER (WATTS): °°,27.DDE*:Power (3) 
ike CHIP uNOGame “St cmenn@ eX ene? 7 POE hie )1( 19). 120) 
"SX, °*POWER (WATTS): °°,22.DDE":Power(4) 
"1X, °°CHIP NOS: °°,5(22.DDE.2X),2X.72.D0E "371249, 7( 25), 1126) 
"SX,°°POWER (WATTS): °°,22.0DE";Power(S) 
Xr "CHIP NOG mee S577 0NE 20.2%. 22-DDE- 71030) .1131).1(32) 
"SX, °*POWER (WATTS): °°,Z7Z2.DDE*:Power(6) 
wie SCHIP NOs meee? aD0EeeX)), 2X, 77.008". 1(36),11 37), teed) 
"SX," "*POWER (WATTS): °°,.22.DDE°:Power(7) 
"tX,°*CHIP HOS: °*,S$(72.DDE.2X).2X.2Z.D0E 727142), 1143), 1144) 
"SX,"*POWER (WATIS): °°,2Z2.DDE°: Power (98) 
"1X, °°CHIP NO9: °°,5(2Z.00E,2X).2X.Z2Z2.DDE°:7T(48),7(49),T(50) 
"SX,°°POWER (WATTS): °°,22.DDE":Power(9) 


"SX, "HEAT EXCHANGERS TEMPERATURES: RIGHT LeFt°** 
SOX se UT TOMS eg alixe 22 ADDENoKe 6 1 (98), 1CS7e 


~wOMee = 1OP3 >" ,24X,.2072.0DE. 2X0" 21160) .1(59) 


"SX,° "BACK PLANE TEMPERATURES ARE :°°* 


MON ess) > «ck, (22.00E. 2X)" 
ml Nats p1S602 > eX. (22. 00E Gene 
Bue 057). > 2X, CLL. DOL Ceuie 


Mle 3). eke (22 -D0E.2K) = 1172) 
ml nemelt 7 4)2°".2%.022.00E.2x)° 51173) 


Slow 75) 62-4 2k.(272.DDE 2X) = 41074) 
ere 07O0>> 62K, (72.006,.2X)° (1075) 


SOA e177) 5°", 2h .0772.0DE .2xX)° 3 TC76) 


ASSIGN @File TO # 
END 
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Calcpiel 


' File NAME $ CalcDivew 

| EDITED BYoyet. 7 BENSDMGi Usk 
| 

! 


THIS PROGRAM ANALYSES THE DATA READ FROM 
' A DATA FILE DES@NATED BYSTHE OPEkeTGr 
! REDUCES THE DATA TO CALCUALTIONS OF NET 
{ POWER, RAYLEIGH NUMBER AND NUSELT NUMBER. 


con iCo/ DCT) 


DIM Emf (70) .Power(9),7(60).Tava(9).Ts(9) 
DIM Tfilm(9) .Qnet(9) ,H(9) ,K(9) .Rho(9) .Cp(9) 
DIM NCS) .Nul9),Ra(9) Delt 3) (ale alon er co 
DIM Gr(3).Betat(3) ,Doowl3) -DtscS) 


! CORRLEATION FACTORS TO CONVERT Emf TO REGREES CELCIUS. 
DATA 0.10086091, 25727.59.-767345. 8.7 euzood5, 
DATA -9247456589,6.98E11,-2.66E13.3.94E14 


READ D(*) 
j 


! RESISTANCE OF THE PRECISION RESISTOR IS: 1% 
Rp=2.00 
i] 


PRINTER IS 701 


INPUT “ENTER THE NAME OF THE FILE CONTAINING DATA" ,Oldfiles 

PRINT USING “15X,""THE RAW emf DATA ARE FROM THE FILE: w", 10A" ;Oldfiles 
INPUT “ENTER THE POWER SETTING “,PowerS 

Soe USING "17X,"" THE POWER SETTING PER CHIP WAS: "=", 10A";PowersS 

PRIN 


BEEP 


BEEP 

ASSIGN @File TO Oldfiles 
ENTER @FilesEmf(*) 

j 


1 CONVERT Emf TO DEGREES CELCIUS 
FOR I=0 TO 60 

Sum=0 

FOR J=0 TO 7 
Sum=SumtD( J) *Emt (1) °J 

NEXT J 

TCI) =Sum 

NEXT I 

i 


d eee! Emf TO POWER 

VoltzEmf (61) 

FOR I=62 TO 70 

Power(J) <Emf (1) *(Volt-Emf(1I))/Rp 
Dpow(J)=.01#Power (J) 

J=J+1 

NEXT I 


1 AREA OF THE BLOCK FACES 
Acen=1.92E-4 


90 


530 Alef=1.44E-4 

640 Aria=!.44E-4 

650 Atoo=4.8E-5S 

660 Abot24.8E-5 

670 Atot=5.76&-4 

630 

690 ICALCULATE THE AVERAGE TEMPERATURES OF THE BLOCK FACES 


710 Tava(t)=(T(0)*Acen+T(1)*#Atop+T(2)*Ariat+T(3) «Alef +T(4)#Abot) /Atot 

<0 Tava(2)=(T(6)*Acen+T(7)*Atoot+T(8)*#Arig+T(9)*Alef+T(10)*Abot)/Atot 
730 Tavg(3)=ti(l2d*AcentT(13)*AtoprT (14) *ArigtT (15) #Alef+T(16)*Abot)/Atot 
740 Tavg(4)=(T(18)*Acen+T (19) *Atootli 20) *#ArigtI (21) *#Alef+T(22)*Abot)/Atot 
750 Tava(5)=(7(24) #Acenti (25) *Atoptl(26)«Arigtl (27) #Alef+T(28)*#Abot)/Atot 
760 Tavg(6) =(7T(30) *Acern+T(31)*#Atop+T (32) #ArigtI (33) *#Alef+T(34)*#Abot) /Atot 
770 Tava(7)=(1T(36)#Acen+T (37) #Atop+T(38)*Arig+1 (39) *#Alef+T(40)*Abot) /Atot 
780 Tavg(8)=(T(42)*#Acen+T (43) *#Atop+T( 44) *#Ariar+T(45)*Alef+T(46)*Abot)/Atot 
790 TavatQ)=( 1048) *Acen+ (49) *Atop+ 1 (50)*Arigth (Si) *Alef+1(52)*Abot) /Atot 


a ; ce est OF THE THERMOCOUPLE READINGS 
cos ieee 

830 Heys 1 

840 

350 RESISTANCE OF PLEXIGLASS AND UNCERTAINTY 
850 ic 764.65 

880 

890 : CHIP BACK SURFACE TEMPERATURES 


900 Ts(1)=T(5S) 
Shit Ts(2Z)=TC11t) 
920 isk =nci 7) 
930 Ts(4)27T(23) 
940 Ts(5)=T(29) 
S50 Ts(6)=T(35) 
960 Ts (737 (41) 
970 1s(8)=1(47) 
980 Ts(9)27T(53) 
990 Tssum=0 

1000 FOR J=1 TO 9 
1010 Tssum=Tssum+Ts(J) 


1020 NEXT J 

1030 ! 

10450 Tsavg=Tssum/9 
105 


1060 ! CGNDUCTION LOSS CALCULATION AND UNCERTAINTY 
1070 Qloss3=(T(17)-T(S6))/Rce 

1080 QlossS=(T(29)-T(SS))/Re 

1090 Qloss7=(T(41)-T(S4))/Re 

1100 Qloss=(Qloss3+Qloss5+Qloss7)/3 

10 Daloss=(Qloss)+*.022 


1 

i 2 ! 

130 ! AVERAGE SINK TEMPERATURE CALCULATION 
140 Tsink=(7T(57)+7T(S8)4+T(S9)+T(60))/4 

1580 Dtsink=(Dt/T(S8)) 

160 L=8.0E-3 

20.) 
180 PRINT USING “1X,"“"CHIP QNET(W) Tava-Ts Nu 
‘UNC IN Nu °",10A" 

130 PRINT 

200 =! 

Za, ! 

220 ! CALCULATION OF NET POWER, Nu, Ra. AND UNCERTAINTIES 


om 


1230 
1240 
1241 
1242 
1244 
1245 
1250 
1260 
1270 
1280 
1290 
1300 
1310 
1320 
l3s0 
1340 
1350 
1360 
1370 
1380 
1390 
1400 
1410 
1420 
1430 
1440 
1450 
1460 
1470 
1480 
1490 
1500 
oie 

520 


l 
FOR J=!t 10 9 
PR NER oto 
ORNs 
PRINT J 
PRINTER [5 701 


1 CALULATION OF Onet 
Qnet(J)=Powerl(J)-Qloss 
Danet=(Dpow!(J)°2+Daloss’2)°.5 


! CALCULATION OF Tfalm 
TEilm(J)=(Taval(J) + Tsink)/2 
Dtfilme(Dtsink 2+Dtava’2)°.5 


1 CALCULATION OF A DELTA TEMPERATURE 
Delt(J)=TavglJ)-Tsink 
Ddelt=(Dtavg°2+Dtsink*2)°.5 

§ 


CALCULATION OF reac COEFFICIENT 
H(J)=Onet(J)/(Atcot*#Delt(J) 
Dh= H(J)*((Dqnet/Onet(J))°* 2+(Ddel t/Delt(J)) 2)°.5 


! CALCULATION OF CONDUCTION COEFF ICEENT 
K(J)=(.65-7.894/74E-4eTfilm(J))/10 


| CALCULATION OF DENSITY 
Rho(J)2(1.825-.00246*Tf11m(J)) #1000 
} 


1 CALCULATION OF SPECIFIC HEAT 
Cp(J)=(.24111143. 7037E-4*TFilm(J)) #4180 
| 


1 CALCULATION OF VISCOSITY 
N(J)=1.4074-2. 964E-2*TFilm(J)+3.8018E-4#Tfilm(J) °2-2.7308E-6*TFilm(J) °3+8. 


1 
1679E-9*TFilm(J)°4 


1530 
1540 
1550 
1560 
1570 
1580 
1590 
1600 
1610 
1620 
1630 
1640 
1650 
1660 
1670 
1680 
1690 
1700 
t7a0 
1720 
1730 
1740 
1750 
750 
1770 
1780 


ets E-6 


CALCULATION OF BETA 
Beta(J)>. 00246/(1.825-.00246*TFilm(J)) 


: CALCULATION OF ALPHA 
ALP at) eR) CROC ee 


CALCULATION OF PRANDTL NUMBER AND UNCERTAINTY 
Pr(J)= NC J)/ALf ald) 

t 

! CALCULATION OF NUSSELT NUMBER 
Nu(J)=HO J) *#bL/K(J) 

Dnu=Nu(J)#*(Dh/H( J) ) 

! 


! CALCULATION OF PERCENT UNCERTAINTY IN NUSSELT NUMBER 
Pernu=(Dnu/Nu(J))*100 

' CALCULATION OF GRASHOF NUMBER AND UNCERTAINTY 
Gr(J)=9.81*Betal(J)*(L°3)*Delt(J)/N(J)°2 
Dar-Gr(J)*(Ddeit/Delt(J)) 


1 CALCULATION OF RAYLEIGH NUMBER AND UNCERTAINTY 


Ra(J)=Gr(J)#*#Pr(J)*1.E-6 
Dra=Ra(J)#(Dagr/Gr(J)) -: 
! 


22 


a CALCULATION OF FLUX BASED RAYLEIGH NUMBER 

QU t 

Te70 Rat (J)=' (9.81 +Beta(J)*h 4*Onet(J)I/(K( JI*N( J) «Al fal J)<Atot) )#1.E-7 

1820 Draf=Raf(J)*(Danet/Gnet(J)) 

1830 =6! 

1340 '! PERCENT UNCERTAINTY IN RAYLEIGH NUMBERS 

1850 Perra=(Dra/Ra(J))#«100 

1660 Perrat=(Draf/Rat (J))<«100 

S70 

cae PRINT BSINGeax 01k, 417%. 22.00E.3 "2. 0metts) velt(s) .Null),Pernu 

pao PRINT USING "6X,""TEMPERATURE BASED RAYLEIGH NUMBER * E-6 IS: °".3D.DDE"; 
a 

1920 PRINT USING "6X.""°% UNCERTAINTY IN THE TEMPERATURE BASED RAYLEIGH NUMBER I 

See 6 ehOuULODDE” :Perra 

1930 FRINT USING “6X. ee CUE RCED RAT ERICH NUMBER E-7 15 Secs. SDePDE” {Ratt )) 
194G PRINT USING "6X,""% UNCERTAINTY IN FLUX BASED RAYLEIGH NUMBER IS: PUD. 

BOE” sPerraf 

1950 PRINT 

1960 NEXT J 

1970 ASSIGN ®File TO * : 

1980 END 
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FastScan 


FROGRAM FASTSCAN 


' 
EUITEO BY El. Ves Bee Der 


| THIS PROGRAM SCANS 3 CHANNELS FOR VOLTAGE MEASUREMENTS 


Ipass*5$99 

Pass*0 

N= 

DIM 711(599).V1(2),Y1(S599) 
DIM 12(599) .V2(2),Y2(S93) 
DIM 13(599) ,V3(2) ,¥3(S99) 
CLEAR 709 

CLEAR 722 

1 


BEEP 

PRINTER [S 701 

eels 

INPUT "ENTER THE FIRST FILE NAME: ° ,NewfilelsS 
INPUT "ENTER THE SECOND FILE NAME: °.Newfile2s 
INPUT “ENTER THE THIRD FILE NAME: °,Newfile3$s 
INPUT “ENTER THE VOLTMETER READING: °.V& 

ae USING *1SX,°° THESE RESULTS ARE NOR STORED ON DISK 
PRINT USING °2SX,°°F ILE: °* ,10A° Newfilels 
PRINT 

PRINT USING °25X,°°F ILE: °", 10A* ;Newfile2s 
PRINT 

PRINT USING "25X,°°F ILE: °°, 110A" Newfile3s 
PRINT 

WAIT 1 


rae ale 

OUTPUT 709; °AET" 

WAIT 2 

Beer 

OUTPUT 722;°T4& FI RI PO 20 ISTI SO! ISTN® 
i 


! LOOP NUMBER ONE 
' 


OUTPUT 709: °AFO ALO* 
OUTPUT 709;°AS°® 

BEET 
Timedate!«TIMEDATE 
FOR Jjy=0 TO Ipass 
OUTPUT 7223713" 
ENTER 72e22:V1(#) 
TIi(PassdVICT) 
Pass=Pass?t] 

NeN+] 

NEXT Jy 
Timedate2=TIMEDATE 
Totaltimel=Timedate2-Timedatel 
OUTPUT 722; °ACO° 
Pass=0 


' LOOP NUMBER THO 
! 
OUTPUT 709; °AFI2 ALI2° 


OUTPUT 709; °AS* 
BEEh 


=) 
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“FASTSCAN’ °° ,10A° 





PLT Pa REEP 

27.4 PR J1=0 10 Ipaiss 
275 Orie eect ° 13° 
c76 PRR 726: V2 8) 
ant T2( Pass) <V2(1) 
<8 Pass*Pass+] 

</9 NEXT Jy 

280 OUTUl 7223 *ACI2° 


281 Pass=0 

282 ' 

a ' tOOP NUMBER THREE 
? { 


28S OUTPUT 709:°AF3t AL3I° 
286 OUTPUT 709:*AS° 


207 BEEP 
288 BEEP 
289 BEEP 


29) FOR Jr=0 T0 Ipass 

292 OUIeEUT 722°:"13°" 

293 ENTER 722:V3( 4) 

294 T3(Pass) =#V3(1) 

235 Pass*Pass+] 

2965 NEXT Jy 

2o7 ! 

298 ' END LOOPS 

299 ' 

300 PRINT USING *15X,°°THE TOTAL TIME ELAPSED WAS (SECONDS):°°,2X,(DDD.DDE)°; 
otaltimel 

301 PRINT 

ane co USING "15X,°° THE TOTAL NUMBER OF SCANS WAS : °°*.2X.(D0DD.D.2X)°:N 
304 PRINT USING °15X.°°THE VOLTMETER READING WAS: °°, ,10A,°3V$ 

306 PRINTER IS 1 


307 BEG 
308 ' 
a TRANSFER FIRST SCAN DATA 


311 CREATE BDAT Newfilel$,20 

312 ASSIGN File TO Newfilel$ 

313 OUTPUT @F ile: Tie) 

314 FOR I1*0 T0 Ipass 

320 Tt(1i)=.10086091+25727 .9471(11)-767345. 8471011) °2+78002556eT1(11)°3 
340 NEXT I. 

341 ! 

a TRANSFER SECOND SCAN DATA 

344 CREATE BDAT Newf:1le2$,20 

345 ASSIGN ®File TO Newfile2$ 

346 OUTPUT sFilesT2(4) 

347 FOR I1*0 T0 Ipass 

348 T2(1i)*.10086091+25727.9472(114)-767345.827T2(11)°2*780025564#T2(11)°3 
349 NEXT Ii 

350 ! 

32), TRANSFER THIRD SCAN DATA 

3$3 CREATE BDAT Neuf ile3$,20 

354 ASSIGN @File TO Newfile3s 

355 OUTPUT @File:T3(#) 

356 FOR I1*0 TO Ipass 

357 T3(11)*.10086091+25727.9813(11)-767345.8#13(11)°2*780025564T3(11)°3 
358 NEXT It 

330 STOP 

400 END 


a5 


4. 


PLO 


FILE NAME se eieGiy 


| EDITED BY Lt T J BENEDICT 


THIS PROGRAM ene THE DATA ACQUIRED BY 


| PROGRAM “FASTSCAN” 
RUNES on 7 0S 


BEEP 

Xmin=0 

Xmax=200 

BEEP 

ee "ENTER MINIMUM AND MAXIMUM Y-VALUES”*,Ymin, Ymax 


Xstep=20 

BEEP 

Ysteps.2 

BEEP 

PRINT “IN:SP1;IP 2000.2000.8000,7000;” 
PRINT ~“SG20 100 02 0e sl= 20h. 
Stx=100/(Xmax-Xmin) 
Sfy=100/(Ymax-Yoin) 

PRINT "PU 0.0 PD" 

FOR Xa=Xmin TO Xmax STEP Xstep 
X=(Xa-Xmin) *Sfx 

PRINT “PAgex.- .03 Xie= 

NEXT Xa 

PRINT "PA 100.0;PU;” 

PRINT “PU PA 0.0 PD” 

FOR Ya=Ymin TO Ymax STEP Ystep 
Y=(Ya-Ymin)#*Sfy 

PRINT “PAD ORs: 1° ¥ita 

NEXT Ya 

PRINT "PA 0,100 TL O 2” 

FOR Xa=Xmin TO Xmax STEP Xstep 
X=(Xa-Xmin) #Sfx 

PRINT “PA°:X,",1003 XT” 

NEXT Xa 

PRINT “PA 100,100 PU PA 100,0 PD* 
FOR Ya=Ymin TO Ymax STEP Ystep 
Y=(Ya-Ymin)*Sfy 

PRINT “PD PA 100,",Y,°YT’° 

NEXT Ya 

PRINT "PA 100.100 PU" 

PRINT sor hm. 2 ese e lio we. 

FOR Xa=Xmin TO Xmax STEP Xstep 
X=(Xa- Xmind*Sfx 

PRINT "PA":X,"7,03 

PRINT "CP -2, -1; re ‘ha; 

NEXT Xa 

PRINT “PU PA 0,0” 

FOR Ya=Ymin TO Ymax STEP Ystep 
IF ABS(Ya)<1.E-S THEN Ya=0 
Y=(Ya-Ymin)*#Sfy 

PRINT "PA 0,°3Y,"" 

PRINT “CP -S.-.25;LB";Ya;°" 
NEXT Ya 

BEEP 

Idi=0 oy 

IF Idl=0 THEN 


96 


610 


63 

640 
650 
660 
670 
680 
690 


710 
720 
730 
740 
720 
760 
ab 
780 
730 
800 


820 
830 
840 
850 
860 
870 
880 
830 
900 
910 
920 
930 
940 
950 
960 
S70 
380 
390 
1000 
1010 
1020 
1021 
1030 
1040 
1050 


BEE 

XlabelS="Time (sec)" 

Beer 

YlabelS="Temperature (C)”" 

Bvereeon 125,42 .FU Pa SUe-Tueer .-LEN(XlapelS)/2:°0:LB°:XlabelS:°° 
cs Siisuren Oso, MENG lapelhiets75° Ul 0. )°LB- -Ylabels:*" 
Pedi CPo0..0" 

bee 

INPUT "ENTER THE NAME OF THE DATA FILE",D_ files 

ASSIGN ®File TO D_files 

BEEP 

Md=0 

BeEEe 

Npairs=600 

BEEP 

PRINTER IS 1 

Sym] 

PRINTER IS 705 

PRINT “PU DI” 

IF Sym=1 THEN PRINT “SM.” 

IF Sym=2 THEN PRINT "SM+" 

IF Sym=3 THEN PRINT “SMo” ‘ 

IF Md>1 THEN 

FOR I=1 TO (Md-1) 

ENTER @File;Xa,Ya 

NEXT 

END IF 

FOR Xaz0 T0 199 STEP .3333333 

ENTER @File;Ya 
Ya=.10086091+25727.9*Ya-767345.8#Ya°2+/8002556*Y¥a°3 
X=(Xa-Xmin) *Sfx 

Y=(Ya-Ymin) *Sfy 

IF Sym>3 THEN PRINT "SM" 

IF Sym<4 THEN PRINT "SR 1.4,2.4" 
PRIN] =Pacek. 1. PD” 

[IF SymoGanen PRINT “SRai.2,1.6°. 
IF Sym=4 THEN PRINT "UC2.4,99,0,-8,-4 
IF Sym=#S THEN PRINT "UC3,0,99,-3,-6,- 
IF Sym=6 THEN PRINT °UCO,.5.3,99.3,- 
Tee syner ener iN UCOL-5.3.,99.,-3 
NEXT Xa 

PRINT. "PU" 


BEBE 
ASSIGN 8File TO * 
END 


oF. 


LO 


thik 
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